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Vibration—Rotation Energies of Planar ZXY,. Molecules 


Part I. The Vibrational Modes and Frequencies 


SAMUEL SILVER* AND EARLE S. EBERs? ** 
Mendenhall Laboratory of Physics, Ohio State University, Columbus, Ohio 
(Received June 11, 1942) 


The normal modes and frequencies of planar ZX Y2 molecules are investigated using the most 
general potential function consistent with the symmetry. The equilibrium configuration has the 
symmetry of the group Cgy and consequently the frequencies are all non-degenerate. By 
developing a set of generalized coordinates having the same symmetry properties as the normal 
coordinates, the secular determinant is factored directly into one third-order determinant, one 


linear factor, and one second-order determinant. 





1. INTRODUCTION 


ERTAIN discrepancies exist between the 
interpretations of the infra-red spectra and 
the electronic spectra of formaldehyde.! Pre- 
paratory to a re-investigation of the infra-red 
data we felt it would be desirable to have a 
complete unified treatment of the vibration- 
rotation energies of planar ZX Y2 molecules. In 
the present paper we develop the normal coor- 
dinates and frequencies; in’ Part II we shall 
consider the quantum-mechanical Hamiltonian 
and its eigenvalues and special problems which 
are pertinent to the interpretation of the 
formaldehyde spectrum. 
The vibrational frequencies have been studied 
previously by Lechner? for a valence force model 
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and by Bernstein* who studied the symmetrical 
modes under the same type of force system. In 
the formulation presented here a general poten- 
tial function is used so that, if desired, other 
force models can be investigated. Also we give 
the normal coordinate transformation which is 
needed for subsequent work. 


2. THE SYMMETRY COORDINATES 


We take as our system of axes the principal 
axes of the equilibrium configuration. The YZ 
plane is the plane of the molecule and the Z axis 
coincides with the twofold axis of symmetry. 
The coordinates of the equilibrium positions are 
designated as xx°, yx°, 2x° (the X atom being 
numbered 1, the Z atom 2, the Y atoms 3 and 4). 
The components of the displacement of the kth 
atom from its equilibrium position are xx, yx, 
zx. The coordinate system and the symmetry 
elements of the molecule are shown in Fig. 1, 
and in Table I are listed the coordinates of the 
equilibrium positions. 


$H. J. Bernstein, J. Chem. Phys. 6, 718 (1938). 
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Fic, 1. The equilibrium configuration and its symmetry 
elements. Ro is the Z—X distance; 79, the X-—Y distance; 
Bo, the ZX Y angle. 


The equilibrium configuration has the sym- 
metry of the group C:,. The irreducible repre- 
sentations of this group are given in Table II. 
Applying the group theory analysis of Wigner‘ 
and Wilson® we find that the normal coordinates 
are all non-degenerate and fall into the following 
symmetry classes: three, totally symmetric, 
transforming according to the irreducible repre- 
sentation A,, one transforming according to By, 
and two transforming according to Bs. The 
secular equation can, therefore, be factored into 
one third-order factor, one linear factor, and one 
factor of second degree. Accordingly, we seek 
three symmetry coordinates belonging to class 
Ai, one coordinate belonging to class B;, and 
two coordinates belonging to class By. The single 
coordinate of class B; is the normal coordinate 
of that symmetry type. 

It is convenient to set up first the coordinates 
expressing infinitesimal rigid translations and 
rotations. These are 


1 4 
2? =— mi Eia 
e M} 2 . 
(a=x, ¥, 23 bi2= &, biy=mi3 bi2= 6), (1) 
and 


; = i? 0 
= y ® m3 (Xigkiy — Xiy kia) 
(1° )b cy 





Qua 


0 0 
(a, B, Y=, Y, OF 3; Xiz=XP 3 Xiy=yeP, etc.), (2) 


4E. Wigner, Géttingen Nachrichten, p. 133 (1930). 
5 E. B. Wilson, Jr., Phys. Rev. 45, 706 (1934). 
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where we have introduced mass-reduced coor- 


dinates fia=mXia. M=>> m; and e. are the 
equilibrium moments of inertia. The choice of 
the principal axes as our coordinate system 
insures the orthogonality of the above coor- 
dinates. Setting Q27=0 and Que=0 is equivalent 
to using Eckart’s conditions’ when the vibra- 
tional motion is referred to rotating axes. 

We require that all the symmetry coordinates 
be mutually orthogonal and also be orthogonal 
to the translational and rotational coordinates. 
The totally symmetric coordinates must be 
linear combinations of 1, ¢, (€3+4), and 
(n3— 4). The coordinates of class By are linear 
combinations of 1, 72, (ns+n4), and (f3— £4). 
The single coordinate of class B, is a linear 
function of the & only. Except for the above 
restrictions, the A; and Bs coordinates are 
completely arbitrary. We have chosen them so 
that, when the ZX Y2 model is reduced to the 
symmetrical X Y3 model, they reduce to the 
symmetry coordinates developed for the latter 
model by Silver and Shaffer.’ We give below the 
inverse transformation expressing the &q in 
terms of the symmetry coordinates. To get the 


TABLE I. Coordinates of the equilibrium positions. 











k\. xp yK° 2° 





mzRo+2myro cos Bo 








1 0 0 ee 

2 0 0 2my?o cos Box t Dor) Xe 
3 0 rosin Bo | mzRy—(mx+mz)ro cos Bo 
4 0 —rosin Bo = 








M =mx+mz+2my. 


TABLE II. Irreducible representations of the group C2,. 











E C2 Ov ov 
Ai;2 1 1 1 1 
Bi; % 1 —1 —1 1 
A» 1 1 —1 —1 
Bs; y 1 —1 1 —1 








6 C. Eckart, Phys. Rev. 47, 552 (1935). ; 
7S. Silver and W. H. Shaffer, J. Chem. Phys. 9, 599 
(1941). 
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PLANAR ZXY2z MOLECULES. I 


latter explicitly, simply interchange rows and columns. 


Imy 4 2m 3 
-( _) (2;°—2:) Qe: B= () (2:9 —23) Qe: etn (Tie) (2:9—2,9)Q5. (3a) 


Str 














MI MI° 2MI° 
yu uy 
2mymz\} —_ mxI,, , . 
m= ( ae) (23° — 22 IS (F 21°S5; 
YY I 
2mxmy ali: 
sania A i 2,9 —3.),S °S.: 
N2 (= ) 23°) + (7 “Y's 5 
uy yy a 
3 Pe mz' sin Bo mz+2my } 
f uF —cos Bo S3 
4 [2(mz+2my cos? BP 2(mz+2my cos? Bo) 








mxmz 0) S 4 ~~. 05 3h) 
+ mI 29° — 21°) S4 pe jo Z3'D5. (3b) 
yy 


zr yy 





‘mz+2Qmy]} 
| ~~ “| Se: 











M 
2my} cos Bo 5 | MxM z Js 2(mymz') sin Bo s 
~ [2(mz+2my cos? Bo)? : M(mz+2my) [2(mz+2my)(mz+2my cos? Bo) }} ' 








re mz} Cos Bo ; mxmy 3 
wi [2(mz+2my cos? Bo) } ; ee 


Mz sin Bo } i 
_ sa(t)'s. 00 
[2(mz+2my)(mz+2my cos? Bo) |! 21° 





Si, Ss, Ss are the totally symmetric coordinates; they reduce to Qi, So, and Sy. of the X Y3 case, 
respectively. S, and S; are of class Bz and reduce to Sy, and Sy of the X Y3 problem. Qs is unique; it 
corresponds to Q3 of the X Y; molecule, representing the bending of the plane of the molecule. The 
coordinates are all normalized to unity. It is readily verified that, because the symmetry coordinates 
are orthogonal to the translational and rotational coordinates, the £4. satisfy Eckart’s conditions. 


3. THE POTENTIAL ENERGY 
Let 6r;; be the changes in the Z—X and X — Y; distances, 6g;; be the variations in the Z— Y; and 
Y;— Y; distances. The most general quadratic potential consistent with the symmetry is 
2V=Ky6ri2+Ko(6rist 631) +K3(8q2s+8q21) +K s5qs1+2K sr 12(5ris+6r14) 
+2K 66r12(5g23+ 5q24) +2K 767 126¢34-+ 2K 6r136714+ 2K o( 67 136G23+ 57 145G24) 
+2K 10(6r136g24+ 67145g23) +2K 11( 6713 +6714) 5g34+2K 1269235g21+ 2K 13(5¢23+5q24)6q34. (4) 


This does not include the energy associated with the bending of the plane of the molecule. The same 
difficulty arises in all planar molecules, since for the latter type of deformation the 6r;; and 6g;; are 
zero to the first order of infinitesimals. It is necessary to add a special term for this motion. Since 
Qs is the normal coordinate for this vibrational mode, we can write for the potential energy associated 
with the bending of the plane of the molecule 


Va,= (Ree/2)Q¢?. (5) 


Returning to (4), we express V in terms of the symmetry coordinates by calculating the 6r;; and 
6q;; in terms of them and substituting into (4). The labor is considerably reduced by considering 
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the A; and Bz modes separately. In the totally symmetric modes we have 





COS Yo 1 sin (Bo+7o) 
6ri3=O6ri4, 6G23= 6G24= 6712 COS Yo— i t-te 
cos Bo 2 cos Bo 
For these modes the potential energy reduces to , 


” 2 2 2 
2VA1= A1167 32+ A226713 + Ag36g34+ 2A 1267126713 + 2A 1367126934 + 242367136934, 
where 
Au =K,+2(K3+Ky.) cos” yot4Ke COS Yo, 


Aoe=2(Ke+Ks)+2(K3+Ki2) cos® yo/cos? By —4(Kog+K 0) cos yo/cos Bo, 
A33= K4+(K3+ Ky.) sin? (Bo +7o0)/2 cos Bo +2Ki3 sin (Bo +o) /cos Bo, 


cos? Yo 
Ay=2K5—2(K3+Ki) 





COS Yo 
—2K;, +2(Ky+Ki0) cos yo, 
cos Bo cos Bo 





sin (Bo+70) Cos Yo 4 Poe (Bo+7o) 


Ais=K7+(K3+Ki2) 
cos Bo 


+2K\3 cos Yo» 
COs Bo 


sin (Botvo) COS Yo 


sin (Bo+7o) 
Aos=2Ki1—(K3+Kiz) ——— 


+(Ko+K1) 





—_ 2Ki3 cos 0/COS Bo. 
cos Bo 


cos? Bo 
Calculating 6712, 6713, and 6qg34 in terms of S;, S2, and S3, and substituting into (6), we obtain 
2VAr=hirS P+ ko2S2? +3353? + 2h12S1S2+ 2k 1351S3+ 2h23S2Ss, 
2mz sin? Bo 


A 22 ii 
My 








Ay 


1 >" cos? Bo Mz 


mz+2my cos? Bo mz my 



































; 2mz sin Bo 
—2 cos BoAi2—4 sin Bo cos BoAi3+ Ao} 3 
my 
k ~ {Au+A 2 Bot2A Bo} 
= d 4 cos 4\j2 COS ’ 
22 mx(mz+2my) (4411 22 0 1 0 
A 1 | sin? Bo 2my sin? Bo cos? Bo 2(mz+2my) cos? Bo 
8 mz-+2my cos? Bol mz+2my — ws neti my oo 
4my sin? Bo cos B 
~ Bo "Are +4 sin Bo cos BoAiz+4 sin Bo cos” Badal ; 
mz+2my 
k 8 {-—2 BoAut+ BoA 
= —ZMy COS PoAii tmz COS L929 
* [2mxmymz(mz+2my)(mz+2my cos? Bo) } ; 
+(mz—2my cos? Bo)Ai2+2mz sin BoAi3+2mz sin Bo COS BoAos} ; 
1 - 
ki3= 2my sin Bo cos BoAi;— Mz Sin Bo COs BoAoe 
7 (mz+2my cos? Bo)[mz(mz+2my) } 
2mz(mz+2my : 
_ (ms r) sin Bo cos BoA33— (mz—2my cos? Bo) sin BoAi2 
my 
. mz COS Bol mz+2my(1+sin? Bo) ] 
+2[(mz+2my) cos* Bo—mz sin® Bo JAi3— Aost } 
— Mi Pee 
kos {2my sin BoA11+2my sin Bo cos® BoAce 


(mz+2my)[2mxmy(mz+2my cos? Bo) }} 
+4my sin Bo cos BoAi2+2(mz+2my) cos BoAis+2(mz+2my) cos* BoAss} . 
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In the Bs modes we have 





6ri2=0; 6rig= —4ri4; 6g23= — 624; 6q31=0, 
whence, 


2 Vb, =2(K2—Ks)6ris-+2(Ks—K 12) 8q23+4(Ko— K 10) 6r1s6q23. (10) 
Calculating 6713 and 6q23 in terms of S, and S; and substituting into (10), we obtain 
(6) 2 VB, =k Se+hssSe2+2hisSiSs, (11) 
where 
bic M { mz(z2°)? sin? Bon. _K mx . . 
e4= — ha Ks) +—(21" sin yo)?(K3— Ki2) — 22°22 sin yo sin Bo( Ky — Kio) 


myx Mz 











(7) : 
ks, = nD } cos” By( Ko —_ Ks) +cos? yo(K3— Ki.) —2 cos Bo cos vo(Ky— Kyo) } ; 
my 
MIz2 ; 
- ask ‘) | Mz" sin Bo cos Bo(K2—Ks) +mxz2;° sin yo cos yo( K3— Ki) 
my I’ Mx Mz 
—[mz22° sin Bo cos yo+mxz," sin yo cos Bo |(Ky—Kio)}. (12) 
The Simple Valence Potential 
(8) The simple valence potential has found wide applicability as a starting point in the analysis of 
the vibrational frequencies of organic molecules. Accordingly, we give below the relations between 
the valence force constants and the generalized constants in the potential energy of the planar ZX Y2 
molecule. For the latter molecule, the valence potential has the form 
™ 2 2 2 . “ ‘i . 
2V=kxzoriotkxy(6rist+6ris) +Rzxy (68 +682") +hkyxvde*, (13) 
where the deformation constants kzxy and kyxy have dimensions of ergs/(radians)*. Considering the 
A, and By modes separately, as was done above, (13) is readily transformed into the general forms 
(6) and (10), respectively. In this way, we obtain 
Au _ kxz, 
4 tan? Bo kzxy 
Aee=2key = yxy+ =). 
ro 
1 kzxy 
A33=———| kyxy+ , 
ro” COs? Bo 2 
Ay=Ai3=0, 
2 tan Bo kzxy 
Ao3= — a (tert ) ; (14) 
ro COS Bo 2 
and 
Ro cos Bo— 70)" 
PR Tee Pe IR som Mi ’ 
| Roro sin By | 
. kzxy 
K;—Ky=—— 


Re sin? vo 
kzxy(R _ 

) Ky—K y="! =— re) (15) 
i Rero sin Bo 
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4. THE SECULAR EQUATIONS 
Since the symmetry coordinates are orthogonal and are normalized to unity, the kinetic energy 
expressed in terms of them is a sum of squares: 
T=}3(SP+S82+8;*) +3Q0E+3(SP2+S;). (16) 
The frequencies are the roots of the secular determinants 
|Rij—D6,;| =O; A=4q?v?C?. (17) 
It is directly evident that Q¢ is the normal coordinate of class B,, corresponding to the frequency 


(in wave numbers) 


1 
ve=—V kee. (18) 


The totally symmetric coordinates give rise to a third-order determinantal equation ; the frequencies 
are the roots of the cubic, 
i — (Rii+Rooths3) 2+ (Riskoo +h iiks3 + kooks3 — kis — Ri2— kos) 

+Rikosthookisthksski2— 2Rioki3ko3— Riikook33=0. (19) 


Let Ai, Az, As be the roots and Q;, Qe, Q3 be the respective normal coordinates. The transformation 
from the symmetry coordinates to the normal coordinates is then, 


Sr hi lie lisy Qi 
So =F lor loo los Qe (20) 
S3 I31 132 1337 \Qs 


where 











k13(Ro2— Xx) — Riokes Ro3(Ri1— Ax) — Rivkis | oe) a 
= 5 lx = +> lye , (21) 
D, dD, D,, 


lk 
with 
D,.= { [Ris(Ro2— Ar) —Ryoko3 | +[Roa(Rir — x) — Riokiz P+( (Ri — zx) (Ree eae. Ar) wae ki2}} J 


The By, frequencies are obtained from the roots of the quadratic equation 








dN — (kaa thos)A+Riskss— his =0. (22) 
Denoting these frequencies as v4 and v5, we have 
oer (s**) +e]. (23) 
2m 2 2 
The transformation from S, and S; to the corresponding normal coordinates Q, and Q; is 
—— s=TqsQat]asQs, Ss=lssQatlssQs, (24) 
lia=lss= =~ lis= —lss — (25) 


[(Rsas—Ma)? +R? : ~ [(Ras— da)? +82 
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Vibration-Rotation Energies of Planar ZXY. Molecules 
Part II. The Quantum-Mechanical Hamiltonian and the Energy Values’ 


SAMUEL SILVER 
Mendenhall Laboratory of Physics, The Ohio State University, Columbus, Ohio,* and 
Department of Physics, University of Oklahoma, Norman, Oklahoma 


(Received June 11, 1942) 


Nielsen’s treatment of the vibration-rotation energies of a general polyatomic molecule is 
applied to the planar ZX Y2 molecular model. Due to the low symmetry of the model, the 
vibrational modes are non-degenerate, and, if the frequencies are well separated in their values, 
Nielsen’s general formulation can be applied directly. For this case the elements of the secular 
equation are evaluated to include all the first- and second-order terms in the Hamiltonian. With 
the view of applying the theory to the analysis of the fundamental bands of formaldehyde, the 
effect of an approximate degeneracy between a pair of perpendicular vibrations, of symmetry 
types B; and Bz, respectively, is considered. For this case, the vibration-rotation energies of the 
states giving rise to the fundamental bands are investigated to the first order of approximation. 
It is shown how, by a modification of Nielsen's general treatment, the secular equations can be 
set up with little difficulty to include all the first-order terms in the Hamiltonian. The secular 
equations for values of J=0, 1, 2, 3, 4 are given explicitly. 





HE quantum-mechanical calculation of the vibration-rotation energies and the associated 

spectra of planar ZX Y2 molecules has been carried out previously subject to a number of 
simplifying assumptions.? Recently Nielsen? expanded the Wilson-Howard Hamiltonian for the 
vibration-rotation energies of a general polyatomic molecule and delineated a procedure for evalu- 
ating the energies to the second order of approximation. In the present paper the general treatment 
is applied to the planar ZX Y2 model and the constants appearing in the elements of the secular 
equation are evaluated in terms of the coefficients of the normal coordinate transformation, developed 
in the preceding paper (Part I), and the coefficients of the cubic and quartic anharmonic terms in 
the potential energy. 

1. THE QUANTUM-MECHANICAL HAMILTONIAN}{ 


The quantum-mechanical Hamiltonian for the system, arranged in terms according to orders of 
magnitude, is 




















RS kg 
a ( “ th th wee har(™+02) 5 wy, = 2rcr;. (1a) 
rer vy z 
H th 3 hi 3 Fé Fr, P? , 
geiees 4 2 o1Pa,,— 2 2 jai th, n +e |n/ Wi: 
rz uy zz 
hi s p-P:; OL’ » p-P. 
<j = 
+r I (P,P.+P.P,) yy C4Q/ WK: po + }° + rP }+ U;. (ib) 
yy 22 rz uy zz 
th 16h PZ Pr P?) 1/07 7. $2 
E =o , in ee at pa. a = tl - 
T, 2 2 SudiPa, F> 2 as | east €utoat fh +3 P "p t5 )+ U2, (1c) 
rz vy 22 rr y zz 





' This paper was presented at the meeting of Section B, The American Association for the Advancement of Science, 
held in Dallas, December, 1941. Abstract, Phys. Rev. 61, 383 (1942). 


Sis The work covered by this and the preceding paper was initiated at the Ohio State University while the author was 
in residence there during 1940-1941. 


7H. Hi. Nielsen, Phys. Rev. 38, 1432 (1931); J. Chem. Phys. 5, 818 (1937); Phys. Rev. 55, 289 (1939). 
*H. H. Nielsen, Phys. Rev. 60, 794 (1941). 


t Only the essential results are given here. For the details of the development of the Hamiltonian, the reader is referred 
to reference 3 above. 
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plus terms which make no contribution to the second order of approximation. q; is related to the 
normal coordinate Q, by qi= (wx/h)*Q:; px is the momentum conjugate to g,; Ui; and Uz are the 
cubic and quartic anharmonic terms in the potential energy. From the symmetry properties of the 
normal coordinates (and hence, the g,), it is seen that the most general forms for U; and U2 are 


U1 =hc{ B11191° +B 222G2* + 833393" + B11291°G2 + B11391°9s + B22192"Gi + B22392°Gs + B33193°q1 + B33293°G2 
+B 441941 +B 4429°2 +B 4439293 + B55195°91 + B55295°G2 + 85539593 + 86619691 + Boo2geG2 
+B 663967G3 + 8123919293 +8 145919195 + 8245929495 + Bs459391gs}, (2) 
U2=het > Vowke Ga FY 11129192 + 11139193 + ¥222192°G1 + Y 22239293 + Y 333193 G1 + Y 33329392 
Foy s44sqsqs + -vessadsgety11209192?-+-yu1aagi'gs? + yseag’ge- + yinssgi2gs? +-7116691°G6? 
+ ¥223392"Qa" + ¥204492"G4 +7 225592°Gs" + Y 226692°Ge" + ¥334493°G4° + Y 335593 G5 + Y336693' Ge 
+7 44559495" + 44669476" + Y 55669596" + Y 1233919293" + Y 1244919294 + Y 12559 19295" + Y 1266919296" 
+1322919392? + Y 1344919394 +7 1355919295" + Y 1366919396 + ¥ 2311929391" + ¥ 2344929394" 
+ 2355929398" + ¥ 2366929396 + Y 1145919495 + 224592195 + ¥ 3345939495 + ¥ 664596 949s 
+ 124591929495 +7 134591939495 + 234592939195}. (3) 


Although the coordinates gs and gs cannot appear singly in the first degree, the product term q.q; 
is admissible since both coordinates transform in the same way under the operations of the symmetry 


group. 
The operators pz, py, P. are the components of the vibrational angular momentum and are given by: 


fv \3 vs\3 , \3 vs\? 
rok laf Jom -()anfenf 2) n-(2) on] 
k=1 { V4 VE J Vs Vi 


rf ve\3 vie\? 
fz (~) nbs, (~*) Gea, 
k=l L\ % V6 ; 
. T fv \? ve\3 ' 
p:=L fz (=) avon,-(*) Gi Pa, 
hed L\ v6 Vi , 


The constants appearing in H are defined by 
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2my i 
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" 2mymz 4 . . mM } . 
2) n=2|| cos? 4 sities pat | iia aaa 
+m rot J 9— 9) sin Bol | k=1, 2,3 
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MI, ,21°-+m,] ¢2(22° — 23°) 


ToL? DP (mz-+2my) (mz+2my cos? Bo) ]! 





Laulsi, (k=1, 8 3 ’ j=4, 5) 














(7) ‘ 
; = ] ° cos Bol +| oe J 9) 
I (mz+2my cos* Bo) pornren neta) _—- 
mxmzM ; . 
in Bols., (k=1, 2, 3 
‘Gece cos? =" nn ' 


See =Lax, (k=4, 5). 
The coefficients /,; and the equilibrium coordinates 2,° are given in Part I. 











SAMUEL SILVER 


2. TRANSFORMATION OF THE HAMILTONIAN AND THE ENERGY VALUES 


Shaffer, Nielsen, and Thomas‘ have shown that the perturbation calculation can be simplified to a 
considerable extent by applying a contact transformation to the operator H, H'=THT™—, which 
carries terms of first order in H into second-order terms in H’. Details of the method are to be found 
in the paper by SNT; a table of the operators comprising the transformation function 7, appropriate 
for the coordinates used here, has been published by Silver and Shaffer.*® 

Due to the low symmetry of the system, the vibrational modes of the planar ZX Y2 molecule are 
_all non-degenerate. In this case, it is possible to remove all the terms from H, into H,’. The trans- 
formed Hamiltonian is 


Hi =H); Hi'=0; (8a) 
H,! = —he{ tr222P 2! + Tyyyyl y+ Tsesek s+ Tesyy(P oP +P YP 2) + tyyse(P YP P+P2P,) 
+ teeee(P2P24+P2P 2) +t y2y(PyPst+PPy)*} 


6 . PZ 6 . Pr, 6 F P? 1 p? ?,? p? 
. | (= a ) pe (= si ) pt (= ma) TP | Phe Tp Tp 
wx uy rz yy zz 


zz 











2 2 
3 Bel 30,902 3 6 5 6 
—he DL a ee +54 [Hic | LD veveeGitD DL (vice — Ninedq?qe? | 
k=1 


es j=1k>j 


2 
$8 B(1—5ix) 
hel > a [lane n2)qA-+40ingPa2+2v2472h0,/ we} 


j=1k=1 2 v;.(4 v;? = v4") 





2 
6 Bi; 9 
= ie] > _ [(2nt—mta,'+4emgrae-+209ipai/h*]], (8b) 


j=4 k=1 2,.(4v;?— v;,°) 


1 (" Pra 4a ) 
Trarz— a 
8hcI™ \ ww? we? w3” 
zz 
1 b? Pred Phat ) 
Tyyyy— : ’ 
8hcI™ \w;? we? ws” 
yy 


1 (= +24") 
Tzzz22— ’ 
8hcI™ w1" ws” w3" 
zz 





where 



































1 1b; debe aszbs 
ne Bhe( 19 hae ar at 
1 bic, dele d3C3 
T yyz2— ~ 33 Oe od T =). 
1 C101 Cole C303 
eer Sy ere we? ). 
1 eg es? 
we Bhc(T PP)? “+"), (9a) 


4W. H. Shaffer, H. H. Nielsen, and L. H. Thomas, Phys. Rev. 56, 895 (1939). : 
5S. Silver and W. H. Shaffer, J. Chem. Phys. 9, 599 (1941). In addition to the elements given in this table we need 
a term in S, 


= Ain (wi? —_ w? iit Wm?) Pididen +o (w? = Wm -™ wp) HE +wm (wn? —we— wot) AEP = 20.010mPrPiPm jie / 


2 
(cp? +1 Font — 2epPeoi? — 2wPomn? — 2wn2"’) 








to remove the @¢ingiqigm term from H. 
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— ’ 
he 82°C? v6 j=1 2hhw;} 

















Dick _ Evk 3 (1426 j.)b Bre; i663 11+ v6") 
he 82°C", j=1 2h ty 3 82°C?» ( v2 — Ve") 
Des: Eyk - —— 
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he 827°C? V6 j=1 2hhw;} 82°C" v6 j=1 ver— v;* 
Ckk €zk : (1 +26;x) 
— a CiBurjs (k=1, 2, 3) 
he 8r°c?v, j=1 2hhw ;3 

2 
Chek €2 3 C;B Kx; (3042 + v6? 
Cok _ k 1+ ‘Buri - rs id  (k=4,5) 
he 89c?8y, = ja 2hiw;? = 82°C? 4 (v4? — v6") 

2 

Cea Eb 3 ¢;B66; €j0(3r6?+ v7) 








cate 4 3 it 


‘eae ’ 
he 82°C? v6 j=1 2hhw;} j=4 82°C? v6( v6? — v;”) 


Bissh1¥20s “ Bessvevevs h ; oe h 
go =———_ +> + ae > re 
4D 03 kat 4Dye45 1287°¢ p21 wi” 








SBruaBiie | $BiiiBeui BreiB iit, Bis3vi( v2 — 442 — v;*) 




















, (k=4, 5) 


, (k=1, 2, 3) 


6 
reg _. Skk 5 


167°C pa1 


Dyin _ (A+ yy + Vn! = 2,207 ~ 27 Vn" _ QV m7 v7) ’ 


(k—1, 2, 3; j=1, 2, 3;1=1, 2, 3) 


Mj; = , 
ines Vk Vj VI 2D 123 
Neu 3 (14261 )BiuiBoar , Bessrs(ve— mt — re) (k=1, 2,3) 
i=1 v1 2D 445 
3 (14-2641) BixiB on a 
him ( Kt) Brat ost _, Peas a(v4? — V4 5  (k=1,2,3) 
l=1 VI 2D 45 
3 (14-2641) Bex 
a ( md PrsBeas (k=1, 2, 3) 
l=1 VI 
2 
— > ene — v4’ — v5") | 
l=1 VI 2D 45 
3 
Nus=¥ Bribes: (k=4, 5) 
l=1 Vi 





(9b) 


(9c) 


(9d) 


The general procedure is to set up the matrix of H’ using as basic wave functions the harmonic 
oscillator functions yv;,(q;) corresponding to vibrational quantum numbers V;, and the symmetrical 
top functions ¥r°(J, K, M). H’ is diagonal in the vibrational quantum numbers (to the second 
order of approximation) and in the quantum numbers, J and M, and has only (K|K), (K|K+2), 
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(K|K+4) elements in the quantum numbers K= —J, —J+1, ---, J. As is well known, the con- 
sequence of the latter fact is that the sub-matrix for a given vibrational state and J value can be 
reduced to a pair of sub-matrices corresponding to even and odd values of K, respectively. It was 
shown by Wilson® and in more detail by Shaffer and Nielsen’ that by taking suitable linear com- 
binations of the yz°(J, K, M) each of the latter sub-matrices can be further reduced. The reduction 
indicated by Wilson can be effected also in the present case where (K | K+4) elements are involved. 
Accordingly, we take as our set of basis functions 


¥ex=II vvi(qgi)Ax= relat —K, M)+wWe(J, K, M)}K>0, 


i=l 


v= YVi(Qi )Ao=Wvvr (J, 0, M), (10) 


tx=II vlan = Vol al(T —K, M)—¥al(J, K, M)jK>0. 
i=1 


In terms of these we have 
(W,*| H’| WV.) =(¥.*| H’|¥,) =0 (11) 


and the following non-vanishing elements: 
(V, Ax| V, Ax) 
(V, ®x| V, &x) 
(V, Ao| V, Ae) =(V, Ao| V, Ao) =(Rat+4Rs)[2f(f—2) }}, 

(V, Ax| V, Ax+2) 
(V, ®x| V, ®x+2) 


= (Rot R2K?+ RK") +6x1(Rit+2Rs)ftsx2Ref(f—2), (K=0,1,---, J) 


= (Ret ReLK?-+(K-+2)")(f-K(K-+1)} (f—(K+1)(K+2)} } 
+inRILF-DG-9}, (K>0) 
(V, Ax| V, Ax-2) 


(V, &x|V, eae = (Ri+RsK?+(K —2)*)) [{f-K(K —1)} {f-(K—-1)(K—2)} }! 


+5xsRefL(f—2)(f—6)}*, (K>2) 
(V, Ao| V, As) =(V, Aa| V, Ao) = ReL2f(f—2)(f—6)(f—12) }', 

(V, Ax| V, Ax+s) ) 

(V, ®x| V, x+4) 

where f=J(J+1), and 


Y, 
R= 441) s{it25% pvt | s{it22 (vit) 


=RL{f—K(K+1)}{f—(K+1)(K+2)} {f—(K+42)(K+3)} 
x {f-(K+3)(K+4)}]!, (K+4>5, K-—4>1) 











kai I” ray 
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+— (tTerzet Tyyyy 2rezyy— 2Tyzyz) } —D,J?(J+1)* 
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XxX, Y, 
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6 E. B. Wilson, Jr., J. Chem. Phys. 4, 313 (1936). 
7™W. H. Shaffer and H. H. Nielsen, Phys. Rev. 56, 188 (1939). 
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with 
h hi 
Z.= 2.70 5 Dag = Tyeys t 3 (Tyyezt T2222) — ET 2eyy— §(Trz2r2tTyyyy) } , 
8r°cl 8x‘ 
zz 
4 
R= — | Tyzyz +3 (Tyyzet Tss22) — Tezyy— 27 2222 — i6(Tr2z2t Tyvuu) I 


aa 5h! hi 
Ry= 4 { [1425 i 4 rity) |-s]i423 ie map ty |- 64 rr —( T2223 — Tyww) J (J +1) ; 
TT 








et if fies al 640 A 
h' 2 
R;= 1282 ——{(trrzz— Tyyyy) —2Tezex tT 2T yyze FAT yzyz} ’ 
hd 
Ro= Saat Teew— 3 (Taeeet Tyvey) | - (13) 
The vibrational term in Ro is defined by 
E, 
Oe ys (Vi 4+DntE Gri( Vict 3)? 1¥ > Gii( Vit d)(Vit3), (14) 
' he k=1 k=l k=1 j>k 
where 
2 2 2 2 
i Bia ee Bir2¥2 Biisvs" Boo1V1 Bo23¥3 B331¥1 
a-ge-— ft —+- + “| + + + + : 
— 16 yy VE 3 z ia 16L4y;2—v.? 4y;2?—v32 4y.?—v,? 4ye?—v3? 4y32?—v? 
Baars a. Bisame |- Re . 2 2 Fe Le 2 2 
i at fe5) —— Le See——(Sas tH fs8) ; 
dnt at — pe? +h aa 2 2» Saat Ses) — 2 ¥ tn 2 (Sat $00) 
3 15 Bins : (8 ne —3 vi”) 2 
Ger =—VEeKkk —— —— 1 — 6...) - (221,2,3 
kk , Iekekeke Yes 2 ( kt) bit ( ) 
3 . (8y,2 —3»;*) 2 
Gi =—Vexkk — >| ———————Biz;;  (k=4, 5, 6 
kk , kklek 2 4y(4n2—vph ( ) 
2 2 
Biixv; rk iVk ; 
Ga=rime— Nine — — |: (k, 7=1, 2, 3) 
4v?— yp? 4y2—y/* 
Bis Vie Vi; MK ; 
G j= jjxk — Nj jxk -——— > stent : “i 47) ; (j= 1, 2,3; k=4, 5) 
4y,2—p;? VE Vj 
; ga=Xe; gs=Ve3 ge=Ze. 
G45= 74455 — Naas; 
VE: V6 
Gre=Yieves— Nixes tZe cie(“* +"), (k=4, 5). (15) 
V6 Vi 
The energies (in term value form) are the roots of the secular equations 
ead | (Wx | H’ | Vo, x) —€dxx-| =0 (16a) 
| (Wr x | H’|V., x) — edxx:| =0. (16b) 


It is evident from Eqs. (12) that (16a) and (16b) can each be factored into sub-determinants cor- 
responding to even and odd values of K, respectively. Illustrations of this reduction are given by 
Shaffer and Nielsen.” 








3. TREATMENT OF AN APPROXIMATE 
DEGENERACY 

The planar ZX Y2 model is of particular 
interest in connection with the formaldehyde 
molecule. Ebers and Nielsen® have attributed the 
complex structure of its infra-red absorption 
spectrum in the 7.5—10y region to a Coriolis type 
interaction between a pair of perpendicular 
vibrations,tf one of class B,; and one of class Bz. 
Their point of view is that such interactions may 
become prominent when two fundamental fre- 
quencies happen to be nearly equal. Nielsen*® 
has considered the effect for a restricted model, 


A B 
A |AA 


















































Fic. 1. Matrix of the transformed Hamiltonian when an 
approximate degeneracy exists between a pair of vibra- 
tional frequencies. 


assuming formaldehyde to be a symmetrical top 
and neglecting the anharmonic terms in the 
potential energy. Wilson® has given a more 
general discussion for the asymmetrical top, but 
likewise neglected the anharmonic terms and 
other perturbing influences. He derived certain 
symmetry conditions which must be satisfied by 
the wave functions for the Coriolis perturbation 
to be significant. The states associated with the 

8 E. S. Ebers and H. H. Nielsen, J. Chem. Phys. 5, 822 
(1 Nhe abstract in reference 1 contains an error of state- 
— Se the pair of modes considered in this 


9H. H. Nielsen, J. Chem. Phys. 5, 818 (1937). 
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frequencies noted above satisfy these conditions. 
In the present section we present a more com- 
plete treatment for the lower energy states 
associated with these frequencies, taking into 
account all the first-order terms in the Hamil- 
tonian. We denote the frequencies as v5 (class By) 
and v¢ (class B,), respectively. 

The effect of the accidental degeneracy makes 
itself felt in applying the contact transformation 
to remove from HZ, the {55 term in the p, com- 
ponent of the vibrational angular momentum 
(Eq. (4)). The element in the transformation 
function 7 required to remove this term is® 


P.{ (ws? + we?) 9596+ 2wswepa,ha,/h? 
56 . 


Pe (wswe) * (a6? — ws?) 





When the denominator becomes small, require- 
ments as to orders of magnitude make such a 
term inadmissible. However, the contact trans- 
formation can be applied to remove all other 
terms in /7;, and we obtain in this case, 


Hy’ =H; 


fseP Al /vs\? ve\? - 
H,'= —_ | (=!) abs,—(“") abo} (17) 
V6 V5 


H,' retains the same form as in 8(b) 
except that we delete the {55 terms 
from the definitions of css and C¢e¢. 





(18) 


Using the set of wave functions defined in Eq. 
(10) as basis functions, we proceed as in Section 2 
to set up the matrix of H’. As before, HH’ is 
diagonal in J and VW, ---, V4. The infinite 
sub-matrix for a given J value and set of vibra- 
tional quantum numbers Vi, ---, V4 has the 
form shown in Fig. 1, corresponding to the 
totality of values of Vs; and Vs. The diagonal 
elements A, A’, A’’, --- arise entirely from JJ)’ 
and H,’ and have precisely the same form as the 
sub-matrices of H’ for given vibrational states 
developed in the preceding section. The defini- 
tions of Ro, ---, Rs carry over completely except 
for the modification noted in (18). The non- 
diagonal elements are contributed entirely by ///. 
Elements A’A’” connect the pairs of states 
between which the approximate degeneracy 
exists, i.e., Vs, Ve and Vs;+1, Ve*F1. Elements 
B’, B”, --+ connect the states V5, Vg and V;+1, 
Ve+1. The ratio of the orders of magnitude of 
the B’ clements to the A’A” elements is 
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(vs—ve)/(vst+ve6), while the contribution of the 
elements B’ to the energies of the lower states is 
of order of magnitude (vs;—v¢)?/v*. We are 
interested here in the ground state and the first 
excited states which give rise to the fundamental 
bands vs and »¢. To the first order of approxima- 
tion, we may, therefore, neglect the elements 
B', B’, ---+. Accordingly, the energies of the 
ground state are given by the developments in 
Section 2, setting V;=0. For the first excited 
states we set up the matrix using the set of 
functions 


Voa(K) =,(Vs=1, Ve=0, K) ; 
Vo(K)=W¥,(V;=0, Ve=1, K), 
Wua(K) =¥.u(Vs=1, Ve=0, K) ; 
Vuv(K)=V.(Vs=0, Ve=1, K). 


(19) 


The elements of A’ and A’”’, as noted above, are 
given by Eq. (12). Since 
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the only non-vanishing elements of A’A” and 
A" A’ are: 


(Wa | Ha! | War) = (Waa | Ha! | Woo) = — (Wee | Hy’ | Vga) 


= —(W,»|Hy'| Vie) =iR, Khe; (21) 


where 





seh [ Vet | (21a) 


on 82°c]® L (»5%6) ; 


From these results and the discussion in Section 
2 of the elements entering into A’ and A” it is 
evident that the matrix can be reduced to two 
sub-matrices involving elements of Vj. with Vis 
and W,, with V,., respectively, and that each 
sub-matrix can be further reduced into a set of 
two according to odd and even values of K, re- 
spectively. The secular equations for J/=0, ---, 4 
are given below: 





P.A,= —Kh?; P:Ap=0; vain e=Ru; €=Row. (22) 
(20) J=1 
P= — Khaz, e€=Roa; €=Roo. (23) 
f— [RoatReatRsatRo+Rot+Rsot {Raa Ryt+2(Rsa— Roo) if Je 
+[Roat Rat Rsat(Riat2Rsa)f J[Rot+RotRss—(Rv+2Rs)f J—Re=0, (24) 
&—[Roat Ret RsatRoatRwotRso— { Raa— Rap t+2(Rsa— Rsv) if Je 
+[Roat Reat R3a— (Ria t2Ria)f JL Ro+Ro+Rse+(Riwt2Rss)f J—R2=0. (25) 
ie Equations (24) and (25), and 
Roa— € (Raot+4Roa)[2f(f— 2) }} 0 
(Riat4Roa)[2f(f—2) ]}! Roat4Reoat16R3a 2iR, 
+Ref(f—2)—« =0. (26) 
0 —2iR Roo t+4R2+16R3p 
— Rwf(f—2)—e 
Equations (26) with a and bd interchanged. (27) 
J=3 
Equations (26) and (27), and 
Root Reat Rea [Riat10Rsat Roof | aR, 0 
| +(Rut2Redf—e — X{(f-2)(F-6)}}! 
[Riot10RsatReaf | Roat9Rea 0 3iR, 
x{(f-2)(f-6)} +81Rs.— —" 
—1tR, 0 RotRot Rs [Rwt+10Rs— Roof | 
—(Rw+2Rs)f—e xX {(f—2)(f—6)}! 
0 —3iR, [Raet+10Rs—Reof ] Rov t9Rw 
X {(f-2)(f-6)}} +81R3,— 
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Equation (28) with a and 6 interchanged. (29) 
J=4 
Equations (28) and (29), and 
Roa—€ (Ria t+2Rsa) Rea{2f(f—2) 0 0 
X {2f(f—2)}4 X (f-—6)(f—12) }3 
ma Roat4Reat16R3a (Raa t20R5a) 2iR, 0 
xy 2)}3 +Ref(f-—2)-e  X{(f—6)(f—12)}) 
Real 2f(f—2) (Ria +20Rsa) Roat16Rea 0 4iR =(). 
X (f—6)(f—12)}* x {(f—6)(f—12)}} +256Rsa—e 
0 —2iR, . 0 Rut4Rx+16Rw (Ra t+20Rs) 
—Ref(f-—2)-e  X{(f-—6)(f—12)}! 
0 0 —4iR, (Ry +20Rs) Rowt16R», 
; x {(f-—6)(f—12)}3 +256Rs»—e| (30) 
Equation (30) with a and bd interchanged and R; replaced by — Ri. (31) 


Ria is the value of the constant R; for the state V;=1, V,=0, and the given J-value; Ri is the 
corresponding constant for the state V;=0, Vs=1. The secular equations for the state V;= Vs=0 
can be obtained from the above by setting R:=0 in each case and deleting the subscripts a and b. 
In the latter case only one of a pair of equations (e.g., (26) and (27)) need be used. The reduction 
of the secular equation, discussed in Section 2, effected by the use of the basis functions (10) is 
directly evident in the above equations. 

The treatment above has given more than is required for a first- order approximation since the 
definitions of constants R; contain terms contributed by H2. These can be dropped from the final 
results with little difficulty. It should be noted that when the degeneracy is complete, i.e., »5— vs=0, 
the non-diagonal elements B’, B’’, --- vanish and the secular Eqs. (22)—(31) give the energies to the 
second order of approximation. 

The application of the results obtained in this paper to the analysis of the formaldehyde spectrum 
is in progress and will be reported in a later publication. 
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The Near Ultraviolet Absorption Spectrum of Monodeuterobenzene 


C. A. Beck, Department of Physics, The Catholic University of America, Washington, D. C. 


AND 


H. Sponer, Department of Physics, Duke University, Durham, North Carolina 


The analysis of the ultraviolet absorption spectrum of monodeuterobenzene has been given. 
The spectrum is, as in benzene, a forbidden transition which is made allowed by the 603 
vibration corresponding to the 606 «+, benzene vibration. The 0,0 band is calculated to be 


38,124, and it occurs with extreme weakness. 





INTRODUCTION 


HE near ultraviolet absorption spectrum of 

monodeuterobenzene lies in the same 
region as that of benzene.'? This 2300—2650A 
region has been measured and the analysis of 
the bands made with the notation of light 
benzene.* 


EXPERIMENTAL 


The monodeuterobenzene was a sample pre- 
pared by Dr. F. A. Miller and Dr. R. C. Lord, 
Jr., of The Johns Hopkins University. It was 
introduced into the absorption cell in air and 
then the cell was pumped out, with two liquid- 
air traps arranged so as to avoid impurities from 
the pump. The spectra showed the presence of 
some light benzene, the percentage of which was 
roughly estimated from intensities to be about 
15 percent. The light benzene lines can be seen 
clearly in Fig. 1, which represents the first weak 


1W. F. Radle and C. A. Beck, J. Chem. Phys. 8, 507 
(1940) (light benzene). 

*H. Sponer, J. Chem. Phys. 8, 705 (1940) (heavy 
benzene). ’ 

3H. Sponer, G. Nordheim, A. L. Sklar, and E. Teller, 
J. Chem. Phys. 7, 207 (1939). 





Fic. 1. Comparison of first weak 


group in light and monodeutero- 
benzene. 


group of both spectra. Figure 2 shows the gradual 
appearance of the monodeuterobenzene spectrum 
with pressure. All spectra were taken with the 
side arm of the absorption cell in a temperature 
bath, varying from 25°C to —70°C, while the 
main part of this quartz absorption cell of 50-cm 
length was at room temperature. A 5-meter 
concave aluminized grating, in the first order, 
with a dispersion of about 3.4A per mm was used. 
The same all-quartz discharge tube was used 
under the same conditions as in previous work.! 


RESULTS 


Table I contains the frequencies of the 
measured band-edges together with the inten- 
sities and assignments. The intensities are visual 
estimates from a plate with vapor pressure vary- 
ing (in nine equal steps between the above tem- 
peratures), the number denoting that exposure 
at which the band is first visible. Broadness of 
bands, edge variance with vapor density, and 
differing microscope magnification contribute to 
the error in measurement. The main bands, given 
to tenths of cm~! and measured under high mag- 
nification, are probably correct within one or two 
cm~!. Other bands are somewhat less accurate. 
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Fic. 2. Absorption spectrum of monodeuterobenzene at different pressures. 


An effort was made to remove all light benzene 
lines from the table but some weak lines may 
have remained unnoticed. The table does not 
represent a complete list of all the other bands 
measured in monodeuterobenzene. Some 60 weak 
and very weak bands which were measured still 
less accurately are not included in Table I. 

The spectrum resembles very much that of 
light benzene, as should be expected, but the 
appearance of the spectrograms is somewhat 
complicated by the simultaneous presence of 
C;H;D and CeHe. 


ANALYSIS 


The analysis of the monodeuterobenzene 
spectrum should be particularly interesting. 
Although the substitution of a light by a heavy 
hydrogen atom reduces formally the symmetry 
Dg, to Co,, the electronic configuration of the 
benzene molecule remains unchanged. We would, 
therefore, expect the characteristic features of 
the benzene spectrum to show up unaltered in 
the monodeuterated compound. These features 
include the appearance of bands representing the 
0-1 and 1-0 transitions (A% and B%, respec- 
tively) of the 606 vibration which made the 
benzene spectrum allowed, and superimposed 
progressions of totally symmetrical vibrations 
and v-v transitions of non-totally symmetrical 
vibrations. The 0,0 band should not occur in a 
forbidden transition. However, since we do have 
the substitution, we actually have a very slight 
deviation from the benzene case which will 
result in a very slight perturbation. Although 
extremely small, it may be just sufficient to 
bring out the 0,0 band. We should then expect 
superimposed very faintly on this 0,0 band the 
features of the allowed C2, transition. 

The pure vibrational spectra of benzene and 
monodeuterobenzene should, contrary to the 
ultraviolet absorption spectra, be different from 





each other, one representing the selection rules 
for Deg, symmetry and the other one for C2, 
symmetry because here the different masses are 
essential to the symmetry type. This has actually 
been found from Raman‘ and infra-red® inves- 
tigations. Turning back to the ultraviolet ab- 
sorption we then have the situation that vibra- 
tions resulting from C2, symmetry will occur in 
the electronic spectrum which is caused by the 
molecule’s electronic behavior almost completely 
according to De, symmetry. For example, the 
totally symmetrical carbon vibration of 981 
cm! (992 in benzene) should represent all 
strong progressions. Other vibrations which are 
totally symmetrical only for C2, and not for Dg, 
symmetry may occur, although with much less 
intensity. We expect especially carbon vibra- 
tions. 

In analyzing the spectrum we then have a 
twofold task: to interpret the main part of the 
spectrum as coming from a forbidden transition 
A,,—B2, as in benzene with a vibrational struc- 
ture analogous to that in benzene, but with 
additional probably fainter bands resulting from 
the vibrational C2, symmetry of the monodeu- 
terobenzene. Secondly, we must search for the 
possible appearance of a weak 0,0 transition, on 
which may be superimposed very faint bands 
with the characteristics of the allowed C2, 
electronic transition. Finally, a comparison of 
this spectrum with the spectra of the monohalo- 
genated benzenes should prove very interesting. 

Proceeding now with the actual analysis we 
first note that the monodeuterobenzene has a 
great number of weak bands, possibly more so 
than the spectra of light and heavy benzene. 


4 A. Langseth and R. C. Lord, Danske Vidensk. Selskab. 
Math. Fys. Meddelel. 16, No. 6 (1938) ; K. W. F. Kohlrausch, 
Phys. Zeits. 37, 58 (1936); C. K. Ingold and collaborators, 
J. Chem. Soc. (1936), p. 925, p. 971. 

5C. R. Bailey, J. B. Hale, C. K. Ingold, and J. N. 
Thompson, J. Chem. Soc. (1936), p. 931, p. 971. 
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This is in accordance with our theoretical expec- 
tations. 

In analogy to benzene we interpret the band 
at 38,640.8 cm~ as A%. This band appears to 
consist mainly of one very strong component and 
the short edge of this was measured for rather low 
vapor pressure. Here, as in all cases, an effort 
was made to have a narrow region of light appear 
between the cross-hair and the short wave-length 
edge of the bands, which degrade to the red. The 
B®, band probably consists of two components, 
and we have chosen the very strong component 
at 37,520.8 cm. The component of longer wave- 
length which is very weak is not included in the 
table. A band at 37,529 could possibly belong to 
B%, although A% has no corresponding band. 

Just as in benzene, the difference between A°o 
and B%, namely, 1120 cm~', is the sum of the 
value of the same vibration in the ground and 
upper electronic states, which vibration makes 
the transition allowed. Similarly, the bands at 
38,554.7 and 37,433.5 were interpreted as C% 
and D°, according to 1-1 transitions of the same 
vibration superimposed on the A% and B% 
bands. The C%s appears to be single while the Do 
appears somewhat like B%. Choosing the strong 
component of D% (the very weak one is not in- 
cluded in the table) for the calculation of the drop 
in the discussed frequency, we obtain 87 cm. 
Comparing this with the A%) to C% distance, we 
obtain an average of 86.5. This together with the 
1120 as the sum of the frequencies yields 603.2 
cm for the lower state and 516.7 cm~ for the 
upper state of this particular vibration. We will 
use 603 and 517 as values, which gives 86 cm~! 
for their difference. The 603 cm~! compares well 
with the Raman‘ effect whete 602 cm~ has been 
observed. A difference of 88 cm~ occurs in light 
benzene,® while it decreases to 82 cm—! in heavy 
benzene. The 0,0 band can then be computed as 
38,124. A pair of very faint bands appears at 
38,123/38,124.5, which we are inclined to identify 
with the 0,0 band. We assume that the occurrence 
of this band is caused by the very slight deviation 
from the Dg, symmetry, in agreement with the 
theoretical considerations. No such band was 
observed in heavy benzene and its occurrence was 





° A difference of 86 cm™ has been reported in reference 
3, but later measurements suggest 88 cm~! as best value; 
see also the last paragraph of this article. 


very doubtful in light benzene. We estimate very 
roughly that the 0,0 band is about 10* times 
weaker than the A. 

The 603 frequency which makes the transition 
allowed represents in benzene an e+, vibration 
(value 606 cm~'). In monodeuterobenzene, this 
splits into a 8; and a vibration, the first one 
being non-totally symmetrical, and the second 
one being totally symmetrical. As follows from 
the motion during this vibration, the 8; keeps 
about the same value as the 606 vibration in 
benzene, whereas the a; is diminished in value. 
This effect is very pronounced if the substituent 
is of a heavy mass as in the halogenated benzenes 
(e.g., in CsH;Cl the 6; is 611 and the a, is 416). 
It seems that in the A% and B% we observe 
mainly the non-totally symmetrical component. 
The totally symmetrical a; component is known 
from the Raman effect to be 598.5. It could occur 
according to the vibrational symmetries, also, 
although very faintly, superimposed on the 0,0 
band. Because of the closeness of the values of 
the a; and 1, the observance of the a; might be 
difficult, but a comparison between the corre- 
sponding bands of the three compounds shows 
that the bands in all three cases are so very much 
alike, that one may conclude from this that the 
a, component does not show up with appreciable 
intensity. 

The totally symmetrical carbon vibration in 
the upper state which was 923 in light benzene 
and 879 in heavy benzene causes in monodeu- 
terobenzene progressions with separations of 920 
cm~!, The first member is at 39,561 of which the 
violet edge was taken for the measurement. Cor- 
responding progressions appear superimposed on 
B%, C%, and D%» as can be seen from Table I. 
The vibration appears superimposed on these as 
well as on bands not yet mentioned. Since we 
interpreted the occurrence of the very weak 0,0 
band as being caused by a very slight perturba- 
tion of the Dg, symmetry, we might expect bands 
of about the same very weak intensity resulting 
from excitations of the 920 vibration super- 
imposed on the 0,0 vibration. There is indeed an 
indication of such a band at the right position. 
Of course there is always a possibility of explain- 
ing very weak bands by high Boltzmann factors, 
so that it is safer not to give any definite explana- 
tions of such bands. 
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TABLE I. 
—ooOoOoOoOoOoooo ———— —<—— 
6 421 + A’, =A'%+4920 6 39241 A*, = A%+920 
42108 6 39228 
7L? = 42087 6 39221 
6 41770 7 39116.7 G%=A%+4+2X 240 
41427 6 39090.6 
6br 41408 A%=A%+4+3X920 Fo, ? w 39077 0,0+953 
7br 41393 5 39068 
6 41205 6 39059 
7 41196 A = A%+2555 5 39022 
6L? 41164 4 38963 
6br 41111 C= C% +2555 6 33088) G'y=G%— 160 
6br 41083 6 38938 
6br 41040 Aly = A'— 160 4 38913 
6 pnd B", = B'% +4920 4 38891 
6 40986 Abr 38869 
6 40949 Cg = C'%— 160 4 38799 G%y)=G%—2X 160 
6 40918 4 38733 
5 40836 B", =B",+-920 4 38688 
6 40729 10 38640.8 A%=0,0+517 
6 40718 7 38598 
6 40597 6 38583 
7 40520.6 8 38554.7 C% = A%— 86 
7 40506.6 7 38504.1 K% = A%— 136 
8 40489 .6 E°, = E%+920 5 38496 
8 meres) A%=A%+2x920 8 38487.8 
8 40472.5 8 38480.1 A'y=A%—160 
br 40390 C% = C%4+2 920 ? 7 38443.0 B°, = B%+920 
6L? 40373 6 38433.7 
6 40340 £1 =E'o+920? K%=K%Q+2x920 | 7 38421.8 
6 40322 7 38413.2 M%=A%—227 
6br 4031 i} A!, = A%— 160 6 38399 
6 40288 Tbr 38395 Clo = C%— 160 
6L? 40252 sup. M%:=M%+2x920 7 38374.5 
5L? 40108 5L? 38352 D°;=D%+920 
7 40077 B'o = B% +2555 6L? 38343 K'9= K%— 160 
6 40039 G,;=G%+920 7 38335 
5 40027} 7 38326.8 
7 40019 7 38320.8 A*)=A%+2X 160 
6L? 40011 7 38300.8 L°,=L%+912? 0%=A%— 338 
5 39964 6L? 38284 
5br 39917 B"y = B',—160 5 38275 B', = B1,+-920 
5 39894 5 38256 My = M%— 160 
4 39878 G1 =G6%,+920 Sbr 38235 C%)= C%—2X 160 
5 39867 5 38227 
7 39808.0 5 38215 N°, = N%+920 
3 39788 5 38210 
6 39743 5 38196 
5 39675 5 38183 D';=D,—160; K%, ? 
Sbr 39650 5 38174 
7 39622 5 38165 
7 39595.7 A%+955 ?; 0,0+1472 w 38161 A%,=A%—3X 160 
39592 ° 4 38156 
8 39581.3 E%=0,0+ 1458 4 38142 Oo = 0% — 160 
8 39571 0,0+ 1448 4 38130 
9 39561.1 A%=A%+920 vew. 38125 0.0 
9 39552.7 38123 ’ 
7 39473 C%,=C%+920 38118 B*, = B%4+920 
7 39466 5 38114 
5 39437 39596 — 160 3 38108 P%, = P%+4-920 
8 39419.6 E, = E%—160 Sbr 38059 
7 39401.4 A!,;=A,—160 4 38036 
7 39392.6 2 38020 
7 39364 B%, = B% +2920 ? 3 38013 
6 39343 4 38006 
39331 sup. M%,=M%+4920 4 37997 A‘j=A%—4 160 
(39307) C1 =C%,—160 ? 4 37982 0% =0%—2X 160 
6L? 39288 2 37975 
6L? 39271 D°,=D%+2x920 4 37971 
6L? 39263 E*y= E%—2X 160 ? 
6L 
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TABLE I.—Continued. 








37969) 
37961 | 
37948 
37914 
37896 
37884 
37860 
37851 
37844 
37828 
37820 
37801 
37796 
37788 
37767 
37739 
37700 
37677 
37668 
37659 
37636 
37630 
37604 
37578 
37564 
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5 37529 

7 37520 B%=0,0—603 
2 37454 

2 37442 

6 37433.5 D%= B%— 87 
4 aoa L%= B%— 134 
1 37381 

1 37375 

6 S700} 

6 37360.7 B'y= B%)— 160 
4 37346 D' =D%\—87 
+ 37328 

4 37295 N= B%— 226 
4 37280 

4 37275 D',= D%— 160 
3 37270 

4 37256 

w 37226 D'p=L%— 160 
4 sai} 

4 37201 B*,)= B%—2X 160 
+ 37185 P%= B%— 338 
es 
w sont 

w 37041 Bi,= B5—3X 160 
lw 37026 P'9= P%— 160 








L means light benzene. 


The 1008 vibration which is a #;, in benzene 
and becomes totally symmetrical in monodeuter- 
obenzene may be found to occur extremely 
weakly, if at all, superimposed on the A% and 
B%, bands. It also, just as the 920 vibration, may 
occur in connection with the 0,0 band. From a 
careful search of the spectrum and a rough 
estimate of the percentual drop of this vibration 
we conclude that bands in the region of 39,590 
may be associated with the superposition of 
this vibration on A%, particularly the one at 
39,595. This would give a frequency value of 955 
cm in the upper state. The band might, how- 
ever, have another explanation or might repre- 
sent the superposition of two different transitions. 
A faint band may be ascribed to a twofold 
excitation of 955 cm. However, no definite 
evidence could be found for the occurrence of 
955 superimposed on the B%, which may be 
expected from the foregoing. A very weak band 
at 39,077 fits the assignment of 0,0+953. 

The 160 progression of light benzene seems to 
occur unaltered in value in monodeuterobenzene. 
The appearance of the bands, however, is dif- 
ferent. While in light benzene the B'y has the 
same appearance as the B% (namely, a strong 
red component and a faint violet companion not 
quite 6 cm~! apart) and the A's has two com- 


ponents of about equal intensity 4 cm~ apart, in 
monodeuterobenzene there appear always two 
bands of about the same intensity, about 8 cm~ 
apart from each other. It is quite possible that 
this doublet character of the 160 series comes 
from a splitting of the corresponding degenerate 
e+, vibration in benzene, although one might 
expect that such a splitting would probably 
produce more diffuse bands on account of the 
overlapping of the components. 

The H series, which in light benzene owes its 
existence to an excitation of an e+, vibration in 
both the ground and excited states, occurs also 
in the spectrum of monodeuterobenzene, as can 
be seen in Table I. So does the G series which is 
an example of a 0-2 transition of this e+, vibra- 
tion superimposed on A°p. 

The bands at 39,571 and 39,581 resemble bands 
in light benzene at 39,561 and 39,569 which have 
been called the first member of an E series arising 
from the excitation of the other e+, vibration of 
1596 cm in the ground state. There appear in 
heavy benzene two bands at 39,693 and 39,706 
which were taken as the first member of a 
similar E series, but it was pointed out that the 
larger distance of 13 cm~ and particularly the 
larger drop of the e+, vibration in heavy benzene 
made the explanation somewhat doubtful. In 
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TABLE II. The positions of the 0,0, A%, and B% transitions 
for C.He, C.H;D, and Ce6Dz. 











CeHe CséH;D CeDe 
0,0 38,090(calc.) 38,123.8 38,292 (calc.) 
A% 38,610 38,641 38,790 
B% 37,482 37,521 37,712 








monodeuterobenzene the bands are 10 cm7! 
apart and the drop would be between those in 
the other two benzenes, but nearer to light 
benzene. The comparison of the E bands in all 
three benzenes makes the proposed explanation 
still seem to be the best, although it is not 
certain. We have bracketed the bands in the 
table but it is quite possible that only one of 
them qualifies for the proposed explanation. If 
we accept it for one or both bands we should 
expect, just as for the 606 vibration, a splitting 
of the 1596 into a totally symmetrical a; and a 
non-totally symmetrical 8, component which 
give rise to a more complicated appearance of the 
E bands in the mono-compound. It is not prob- 
able that such a splitting occurs; there are two 
bands at 39,539 and 39,553 which fall in the A®, 
region and might be taken as the pair belonging 
to the second component, but since we could not 
detect any indication of a splitting in the A and 
B bands, we hesitate to assume it for the E bands. 
We have not found any conclusive evidence for 
the occurrence of the totally symmetrical com- 
ponent of 1594 in the ground state and that of 
the non-totally symmetrical component of 1576 
is very uncertain because it would coincide with 
B+. 

Bands at 38,413/422 and 37,295 can be con- 
sidered analogous to the bands M% and N% in 
heavy benzene. The separation of the M4 and 
N®, from the A% and B% was 202 cm“ in heavy, 
263 cm in light benzene, and is 227 in mono- 
deuterobenzene. The origin of these bands has 
been discussed previously? and it was pointed out 
that the occurrence of bands displaced by 207 
and 263, respectively, from the 0,0 bands may 
indicate a square combination of symmetry et, 
for these bands. In monodeuterobenzene there 
occurs correspondingly a band displaced from the 
0,0 band by 227 cm~ toward longer waves, 
namely at 37,896 cm~. If we take this fact as 
support of the above-mentioned proposition and 
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if, from intensity considerations and in analogy 
to the benzenes, we estimate the Boltmann 
factor of My as being of the same magnitude as 
that of D'), we may conclude that /7 vibrations 
of about 700-800 cm! are involved in these 
bands. The e~, and e*, vibrations were proposed 
previously as possibilities. No further decision, 
even with our additional data, can be given, 
however. The 920 progression and the difference 
of 160 appear superimposed on both the 11% 
and N°, bands. 

In heavy benzene the letters K°) and L° were 
assigned to bands which appear midway between 
C% and A's, and D% and B'o, respectively. Their 
distances from C% and D% were 113 cm~. Light 
benzene showed no analogous bands of com- 
parable intensity. In monodeuterobenzene we 
have found bands displaced from C%) and D% by 
135 cm to the red. Their intensity is definitely 
smaller than that of the corresponding heavy 
benzene bands. It seems obvious that the bands 
represent 1—1 transitions but we are not able to 
assign them to a definite vibration. The super- 
imposed 920 and 160 progressions do not stand 
out clearly because of some overlapping. 

The totally symmetrical HT vibration which 
was found to be 2565 cm! in the upper state of 
C.sH¢s obviously occurs also in monodeutero- 
benzene in the band at 41,196 with a value of 
2555 cm. It also occurs superimposed on B%, 
and (% . It is not possible to identify the totally 
symmetrical D vibration. The band in which it 
could occur in the ground state must be too weak 
for detection and as for the upper state there are 
too many bands which could be explained with 
the help of this vibration. 

Finally, a few words may be said about the 
fine structure of the bands. As in the benzenes 
many bands appear to be double-headed. (The 
doublet character of the 160 series has probably 
a different cause and has been mentioned before. ) 
We have taken here for the calculations the 
strong head. For the A% and C°y bands this is the 
only measured edge, as also for B% and D%. In 
the 920 progression of the A bands it is always the 
violet head if the two adjacent bands (for ex- 
ample, A®%, 39,553 and 39,561) constitute one 


. band. The 920 progression of the B bands looks 


single. In this way our scheme of differences and 
frequency values was obtained. It is very 
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doubtful whether the faint companions on the 
violet side of the By and D%, 37,529 and 37,442, 
respectively, really are parts of these bands. We 
are inclined to believe that they are not, and we 
would like in this connection to take this question 
up again for the light benzene. 

For CgHg it had been tentatively suggested? to 
consider the corresponding strong red and weak 
violet bands of B%, Do, C% as R and Q branches 
of one band. The single nature of A% was ex- 
plained by assuming that the strong absorption 
masked the structure and that the separation of 
6 cm! found in the other bands between the R 
and Q bands, must be subtracted from the 
measured A °y edge to obtain consistent numerical 
relations. The calculations were based on the 
measured value of 38,614.7 cm~! of A%. Since 
then it has been found in heavy benzene that a 
consistent scheme of relations and frequency 
values could be obtained by using always the 
strong edges and the present investigation yields 
the same result for monodeuterobenzene. More- 
over, a remeasuring of the benzene bands 
changed the value of 38,614.7 for A% to 38,610.7 
cm so that the formerly obtained numerical 
relations do not come out consistently any more. 
This means that the difference between the 
values of the 606 e+, vibration in the ground and 
excited states as calculated from the distance 


A®,—B%) and from the known Raman value of 
606 does not agree any more so closely with the 
distances A%—C% or B%—D%. In fact, this 
close agreement is gone whether we subtract 
6 cm~ or not. Measurements of the bands with 
high dispersion* seem also to furnish doubt that 
the previous procedure of subtracting 6cm™ from 
the A% edge when making calculations is valid. 

From the remeasured bands the 86 cm! dif- 
ference becomes rather 88 cm~!. The present 
small discrepancies would disappear if the value 
of the e*, vibration in question were 608 instead 
of 606 cm-!. The Raman frequency of 606 seems, 
however, well established. Another possibility is 
that the A) is still 2 cm too high, although we 
feel this is a remote possibility only. It is, of 
course, always difficult to make sure that the 
edges of the bands are measured at a comparable 
‘appearance pressure.’’ With increasing pressure 
the bands broaden and the edges move a little 
toward the violet. Finally it is also possible that 
a slightly different structure of A%y and By could 
account for the deviations. Only from measure- 
ments with high dispersion may further informa- 
tion be expected. 


*We are deeply indebted to Dr. Mark Fred and Dr. 
Anthony Turkevich for sending us enlargements of the 
A% and B% bands taken under very high dispersion. 
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Entropy of the Monomeric Forms of Formic Acid and Acetic Acid 


J. O. .HALForD 
Department of Chemistry, University of Michigan, Ann Arbor, Michigan 


(Received May 19, 1942) 


A recent measurement of the entropy of gaseous formic acid at its equilibrium vapor pressure 
is combined with vapor density data to give 60.0+0.3 as the entropy of the monomer at 25° and 
one atmosphere. The corresponding dimer entropy is 83.1+0.3 e.u. Pauling’s residual entropy of 
R In 2 per dimer weight for random orientation of hydrogen bonds in the crystal has been in- 
cluded in the above figures in order to allow a reasonable entropy for the torsional motion of the 
hydroxyl group in the monomer. The same correction has been applied to a previously published 
result for the acetic acid monomer to yield 70.1+1.0 as a revised value. It is highly probable 
that only a single significant potential minimum of undetermined breadth and depth occurs in 































the rotational cycle of the hydroxy! group. 





FORMIC ACID 


CCURATE entropy values have recently 
been obtained by Stout and Fisher! for 
liquid and gaseous formic acid at 25° and the 
equilibrium vapor pressure. Their measurements 
include the redetermination of the vapor pressure 
(43.1 mm) and the heat of vaporization (4754+5 
cal. for 46.026 g). It will be shown that their re- 
sults can be combined with Coolidge’s? vapor 
density data to give the entropy of the monomer 
as a hypothetical perfect gas with a probable 
error less than 0.3 e.u. 

The calculation follows the same lines as the 
one recently made by the writer for acetic acid.* 
It is evident from the vapor density measure- 
ments, as Coolidge has stated, that at pressures 
below saturation, in the temperature range from 
10° to 50°, the vapor is adequately described as a 
mixture of perfect gases in equilibrium according 
to the equation 


(HCOOH)2=2HCOOH. (1) 


At higher temperatures and pressures the equilib- 
rium quotient decreases with increasing pressure. 
The equation 


log k= 24.6070 —3754.9/T — 4.7034 log T (2) 


has been derived by the method of zero sums 
from fifty-four vapor density measurements made 
at temperatures from 10° to 80°. It was found 
that by including some data taken above 50°, the 
precision in determining the slope in the lower 

1J. W. Stout and L. H. Fisher, J. Chem. Phys. 9, 163 
(1941). 


2 A. S. Coolidge, J. Am. Chem. Soc. 50, 2173 (1928). 
3 J. O. Halford, J. Chem. Phys. 9, 859 (1941). 
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temperature range could be increased. This in- 
crease occurs at the expense of accuracy in the 
heat capacity term, which fortunately is of no 
importance to the present calculation. The 
method of least squares gives 


log k= 27.685 — 3904.3/T—5.745 log T.  (2’) 


The sum of the squares of deviations, however, is 
about the same for Eq. (2) as for Eq. (2’). For the 
heat of dissociation of the dimer, the above 
equations lead to the expressions 


AH (dis) = 17,182 —9.348 T (3) 
and 
AH (dis) = 17,866 — 11.399 T (3’) 


and to the values 14,395 and 14,460 cal. at 25°. 
Further application of the least squares method 
with AC> arbitrarily set at zero yields 14,272 cal., 
which is below the limit of reasonable interpre- 
tation, since the curvature is obvious from the 
plotted function. By inspection, probable values 
lie well within the limits 14,400+120 cal. 

The probable error of a single determination of 
log k may be taken as (Zv?/n)}=0.0122, where 2 
is the deviation and n is the number of measure- 
ments. This quantity represents the largest rea- 
sonable estimate of the probable deviation of the 
mean from the true value. A more representative 
estimate would be lower, in the direction of 
(Zv*/n?)!, the quantity used for random devia- 
tions in m measurements of a single point. It 
seems safest here to use the larger probable error. 
Conversion to a probable error in the heat of 
dissociation is made by assuming that a plus or 
minus deviation of 0.0122 in log k might occur 
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with equal probability at either 10° or 80°, and 
applying the expression A7= — Rd In k/d(1/T). 
The result is a probable error of +80 cal. In 
subsequent calculations, the rounded figure 
14,400+100 cal. is used for the heat of dis- 
sociation. 

From Eq. (2), R29s.1.= 2.361 and the degree of 
dissociation a= (k/(k+4p))!=0.1162. Individual 
heats of vaporization for the monomer and dimer 
may be obtained from the equations 


AH(vap) 1-— 
Midian 





AH(dis) (4) 


AH,(92 g) =AH(vap) —aAH (dis). (5) 

For the corresponding entropies of vaporization, 

AH(vap) 1-—a 
+ 








AS;(46 g) = - AFi (dis) 
2ap 
+R In ————__ (6) 
760(1+a) 
— AH(vap) oAH(dis) 
Vi 
AS,(92 g) = ip) a is 
4 (l1—a)p 
-— 
—, (7) 


" 760(1-+a)" 
where AH(vap) and A//(dis) apply to 92 g of the 
acid and is the equilibrium vapor pressure. 

Numerical substitution into Eqs. (6) and (7) 
gives 28.48 and 20.11 for the respective entropies 
of vaporization. By combination with the Stout 
and Fisher result, 30.82 e.u. per monomer weight 
of liquid, the entropies of the hypothetical perfect 
gases at 25° and one atmosphere are found to be 

S:=59.30 and S2:=81.75 (prelim.). 


The sum of probable errors in S; introduced by 
the heats of vaporization and dissociation and by 
the pressure and the degree of dissociation is 0.3 
e.u., and the combined probable error is 0.2 e.u. 
The probable error in S, is about the same, in 
spite of the double weight of acid. 
Consideration may now be given to the ques- 
tion of residual entropy in the crystal, as pre- 
dicted by Pauling.‘ For this purpose the entropy 
contributed by torsional motion of the hydroxyl 
group is found by subtracting the sum calculated 
for the remaining degrees of freedom from the 
monomer entropy. All angles between single 
bonds are assumed to have the tetrahedral 
*L. Pauling, J. Am. Chem. Soc. 57, 2680 (1935). 





value and interatomic distances are taken as 
C-—H=1.09,C=O=1.24,C-—O=1.43 and O—H 
=0.97A units. The mean of moment of inertia 
products for models with cis- and trans-orienta- 
tion of the hydroxyl group, differing by only 0.1 
in their contribution to the entropy, has been 
used. Frequencies of the normal vibrations are 
taken from Bonner and Hofstadter and from 
Herman and Williams® as 3570, 2940, 1221, 1190, 
1111, 1075, and 658 cm-. There can be little 
doubt of the correctness of assigning the lowest 
frequency as a fundamental, and only this one 
must be accurate and correctly assigned for the 
present purpose. 

The calculated monomer entropy, exclusive of 
the torsional motion, is 59.25 e.u., to be compared 
with 59.30 based upon the thermal and vapor 
density data. The contribution of the internal 
rotation or vibration is then 0.05+0.3+.S(res), 
where S(res), the entropy at the absolute zero, 
may be zero, or, if Pauling’s prediction is correct, 
Rin 2=0.69. From these figures the larger 
torsional entropy, 0.74+0.3 e.u., appears by far 
the more probable. This corresponds to a har- 
monic frequency range from 380 to 600 cm~', 
whereas the zero residual entropy would corre- 
spond to the range from 650 to 1150 cm~', which 
is less reasonable. The monomer and dimer 
entropies are therefore revised to 


S,;=60.040.3 and S.=83.1+0.3. 


If the hydroxyl group is arbitrarily treated asa 
symmetrical rotator with a moment of inertia of 
1.35 X 10-* it will contribute in free rotation 4.45 
units to the entropy at 25°, leading to a total free 
rotation entropy of 63.7. This treatment omits 
some small subtractions called for by Kassel’s 
determinant.® It is easy to show, however, that 
the refinements of his interesting and ingenious 
procedure are numerically unimportant except 
for molecules of very low molecular weight. 

The above calculations indicate that in formic 
acid an entropy deficiency of 3.7+0.3 units 
results from hindrance of the internal rotation. 


ACETIC ACID 


It is probable that acetic acid also has residual 
entropy at the absolute zero due to random 


5 L. G. Bonner and R. Hofstadter, J. Chem. Phys. 6, 53 
(1938); R. C. Herman and V. Williams, J. Chem. Phys. 8, 
447 (1940). 

®L. S. Kassel, J. Chem. Phys. 4, 276 (1936). 
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orientation of hydrogen bonds. The value re- 
cently published by the writer should therefore 
be increased from 69.4 to 70.1 (+1.0) e.u. The 
decision is based entirely upon the more accurate 
formic acid calculation, since an unknown uncer- 
tainty in the entropy of liquid acetic acid makes 
it necessary to leave a relatively wide allowance 
for error. For this reason the residual entropy 
problem was not considered in connection with 
the original acetic acid calculation. 

In addition, a revision of the estimated free 
rotation entropy is in order. The vibrational 
contribution was taken from the mechanically 
analogous acetone molecule, for which it had been 
assumed that rotation of the methyl groups is 
practically unhindered. A recent simultaneous 
consideration of acetic acid, acetone and iso- 
butylene, in conference with Dr. R. H. Gillette, 
has led to the conclusion that the vibrational 
entropy of acetone and acetic acid should be 
about 1.2 units higher than originally assigned. 
This raises the free rotation value for acetic acid 
from 72.7 to 73.9, and leaves 3.8 units for the 
deficiency below the free rotation entropy. Of 
this quantity, by analogy with acetone, about 0.6 
should be subtracted for the methyl group, 
leaving 3.2+1.0 for the hydroxyl group. 


THE CARBOXYL GROUP 


The estimated probable errors are large enough 
to obscure any difference between the entropy 
deficiencies calculated for the two acids. It is 
permissible to conclude that the deficiencies are 
large and roughly equal, and to add the suspicion 
that a greater hindrance to free rotation may 
occur in formic acid. 

The above revision of the acetic acid entropy 
produces no appreciable change in the calculated 
entropy deficiency. Consequently the tentative 
explanation made in connection with the first 
calculation is still equally applicable. If hindrance 
of rotation was ascribed to a single sinusoidal 
potential valley covering one-third of the cycle, 
the large deficiency could occur with a potential 
barrier of approximately 3000 cal. For acetic acid 
this interpretation could only be called somewhat 
probable because of the rather large probable 
error assigned to the entropy deficiency. 

The more accurate formic acid result permits a 
more objective consideration leading to the same 
type of description. There must be at least one 
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important potential minimum to account for the 
large deficiency below free rotation (3.70.3 e.u.). 
From the shape of the molecule alone as many as 
three more minima are possible, but they are 
readily shown to be relatively unimportant. For 
this purpose it is convenient to consider the 
contribution of the torsional motion (0.75+0.3 
e.u.) as an entropy of mixing — REX; In X; where 
X ; is the fraction of molecules to be found in the 
ith energy state, vibrational or rotational. If, for 
example, there were two states of approximately 
equal probability, the entropy would be R In 2, 
or 1.38. Still larger values would result from the 
occurrence of three comparable states, or from 
the subdivision of either or both of the two states 
just considered. If all the molecules were divided 
between two states, about 85 percent of them 
would be in the more populated one. 

It is probable, then, that the one minimum 
description is at least qualitatively correct, since 
the occurrence of other important minima would 
cause distribution of the molecules over two or 
more comparable states and the resulting entropy 
would be much larger than the calculated value. 
The results, however, do not permit any specific 
conclusions about the height or effective width of 
the potential valley. The barrier will not be 
unusually high if the valley does not cover more 
than about one-third of the cycle. 

A single configuration might be considerably 
stabilized by resonance of the arrangements I 


and IT. 


O—H 6-H H—O 
| | | 
C C C 


TaN I MS FM if 
H O H O H O H O 
I II III IV 


Formulas II and IV would produce potential 
minima corresponding to cts-trans-isomerism 
which would presumably favor II because of the 
attraction of the negatively charged oxygen atom 
for the proton. If, in addition, a repulsion be- 
tween the hydrogen atoms occurred, IV would be 
still less stable relative to II, and the trans- 
minimum might be substantially eliminated. The 
configuration described in terms of the cis-forms I 
and II would then predominate over the alter- 
native situation embraced by formulas III 


and IV. 
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The mole fraction scale gives simpler (i.e., more nearly 
linear) relations between the variables determining kinetics 
on changing solvents for the bromacetate-thiosulphate 
reaction than does the usual volume molar scale. The 
generality and basis of this finding is discussed from the 
viewpoint of Rabinowitch’s cage model theory as well as 
from that of equilibrium theory. The quenching of fluores- 


1. INTRODUCTION 


ELECTION of the proper concentration scale 

to be used in comparing either reaction rate 

or equilibrium constants becomes a problem of 

importance when the molecular volumes of the 

solvents differ considerably. For example two 

bimolecular rate constants which are equal on the 

mole fraction scale in methyl and in iso-amyl 

alcohol stand in a ratio of 1 to 3 when compared 
on a volume molar basis. 

In calculating thermodynamic properties, it 
has been customary to employ the mole fraction 
scale, whereas in reaction kinetics the universal 
use of the volume molar scale has not been 
questioned until recently.” 

It has been suggested that apparent anomalies 
in the effect of added sucrose and ethanol upon 
the rate in aqueous solution of the reaction! *~7 


S,0;-+ BrCH2COO- = Br~+CH2S:0;COO- 


can be eliminated by the conversion of the data 
from volume molar to mole fraction units. 

In order to investigate this suggestion and to 
attempt to correlate residual and unexplained 


* Publication assisted by the Ernest Kempton Adams 
Fund for Physical Research of Columbia University. 

Constructed from a dissertation submitted in partial ful- 
fillment of the requirements for the degree of Doctor of 
Philosophy in the Faculty of Pure Science of Columbia 
University by Hubert G. Davis in October, 1941. 

1G. Scatchard, J. Chem. Phys. 7, 657 (1939), discussion 
by LaMer, p. 662; Ann. N. Y. ‘Acad. Sci. 39, 341 (1940). 

7. 2s Bell, Ann. Repts. on Progress of ‘Chemistry 36, 
82 (1939). 

3LaMer and Kamner, J. Am. Chem. Soc. 57, 2662 
(1935); 57, 2669 (1935). 

‘Straup and Cohn, J. Am. Chem. Soc. 57, 1794 a 
a — Dissertation, Columbia University, N. Y. 

®°V. K. LaMer, J. Frank. Inst. 225, 709 (1938). 

7 Kappana, J. Ind. Chem. Soc. 6, 419 (1929). 


cence, which involves diffusion as a rate determining step, 
is a more nearly linear function of the fluidity when the data 
are expressed in mole fraction than in the volume molar 
scale. Experiments are in progress to determine, if possible, 
whether the cage model or a modified Smoluchowski diffu- 
sion treatment fits the facts better for low activation energy 
reactions like the quenching of fluorescence in solution. 


solvent effects, we have further studied the effect 
of solvent and of temperature upon this reaction. 
Data on the solubility of thallous chloride in 
solvents of differing molecular volume have been 
obtained for purposes of comparison. 


Ions in Solution 
Electrostatic theory, applied to a reaction in 
dilute solution between two ionic reactants A and 
B of charge Z that form a double-sphere critical 
complex, leads to the well-known rate equation*® ® 


ZaZperNy 1 
DRT rat+re 
ZaLpeNn K 
+ 7 ’ 
DRT 1+ka 


where a is the Debye-Hiickel first approxi- 
mation parameter, ra+rz is the separation 
of the centers of charge in the critical com- 
plex, «= (8re?Nx?2u/DRT1000)! from the Debye- 
Hiickel theory,’® and the other symbols have 
their usual significance. Analogous equations 
have been developed from the Debye-Hiickel 
theory for the constants of ionic equilibria. For 
the solubility of a strong binary electrolyte AB 
the expression is®* ! 


ZtNw {11 
In S=In Sp-»— (-+-) 
2DRT\ra re 


Ze Nw K 
+ ——. 
2DRT 1+xa 


In R=In kp-x — 








(1) 





(2) 


8 Christiansen, Zeits. f. physik. Chemie 113, 35 (1924). 


*G. Scatchard, J. Am. Chem. Soc. 52, 52 (1930); Chem. 
Rev. 10, 229 (1932). 

10 Debye and Hiickel, Physik. Zeits. 24, 185 (1923). 

1 Ricci and Davis, J. Am. Chem. Soc. 62, 407 (1940). 
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In such equations the logarithm of the con- 
stant, or the change in standard free energy of the 
equilibrium between reactants and critical com- 
plex, is expressed as the sum of three terms,® the 
first equal to the value at infinite dielectric con- 
stant and thus free of all electrostatic effects; the 
second, or self-potential term arising from the 
work of charging the ions in a medium of die- 
lectric constant D; and the third expressing the 
effect of the ionic atmosphere at an ionic 
strength yp. 

The third term of Eqs. (1) and (2) offers a 
simple method of extrapolation to infinite dilu- 
tion (kx =0), by a suitable choice of the parameter 
a. When the simple Debye-Hiickel theory cannot 
be applied, this term may be modified by use of 
the idea of ion-association,’ 8 but the extrapo- 
lation is often unsatisfactory. 

Both Eqs. (1) and (2) predict that the plot of 
log k,-o or log S,~o vs. 1/D will be linear. The 
value of r4+7rze or 1/r4+1/rz can be determined 
from the slope. Such a linear plot is often found, 
but there are serious objections to the literal 
interpretation of the self-potential term: the 
value of r4+7z is often either too large or too 
small for the physical model" !* and may even be 
negative; it is in general dependent on both 
temperature and solvent; the use of the macro- 
scopic dielectric constant has long been con- 
sidered questionable, since the orientation of the 
solvent dipoles in the field of the ions will make 
the effective dielectric constant lower.'® The plot 
of log k.~o vs. 1/D is not always linear and indeed 
always shows curvature at sufficiently low die- 
iectrics. Imperfect extrapolation to zero ionic 
strength is an important cause of such curvature, 
but there is reason to suspect that even when the 
extrapolation is valid, a linear plot may not 
always be obtained. 


The Normalized Mole Fraction Rate Constant 


Conversion of a second-order reaction rate 
constant k from volume-molar units, liters per 
mole, to mole fraction units, total moles per mole, 
requires multiplication by the factor, total moles 
per liter. If the mole fraction constant is divided 

2 Bjerrum, K. Danske Vidensk. Selskab, 7, No. 9 (1926). 

18 Hogge and Garrett, J. Am. Chem. Soc. 63, 1089 (1941). 

4 Amis and LaMer, J. Am. Chem. Soc. 61, 905 (1939). 


1 Amis, J. Am. Chem. Soc. 63, 1606 (1941). 
16H. Frank, J. Am. Chem. Soc. 63, 1789 (1941). 


by the total number of moles per liter in a 
standard solvent, a normalized constant, R¥, is 
obtained. This is proportional to the mole frac- 
tion constant, but becomes equal to the volume 
molar constant when the standard solvent is 


used.! 
et > Moles/liter (3) 
—— (Moles/ liter) standara ‘ 


Water at 4°, with 55.5 moles per liter, has been 
chosen as the standard solvent to which the data 
of this paper on the bromacetate-thiosulfate 
reaction have been referred. 

Consider a reaction of the type 





ki ks 

A +Bax — Products, 
where X denotes the critical complex. Let Kx be 
the equilibrium constant (Nx/NaNz)(yx/yays) 
of the former equilibrium, NV = mole fraction and 
y=mole fraction activity coefficient, and let the 
rate determining step of the over-all reaction be 
the unimolecular decomposition of the critical 
complex: 








dj X | NaNp YaYb 
anne X abl: [xX]. (4) 
dt Nx ‘x 
—d[X ]/dt 


k= Molar rate constant = 


[A J[B] 
vars Na Ne [X] 
we 4 [A ] [B] Nx 
=hKr (x =) (5) 


YX liter 











= K3IAX 





where m=2 is the total number of reactants. 
Generalizing the definition of K” (Eq. (3)), we 


have 
moles\ “—-» 
(x i ) 
liter solvent 
moles \ &"~? 
==) 
liter standard 
YAYB 


=k;Kx . (6) 
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Since ks has the dimensions of frequency, we 
denote it by the symbol v. If K, is replaced by its 
thermodynamic equivalent, Eq. (5) becomes for a 


bimolecular reaction 
moles \ ~! 
r ) 7 


liter 





k= pera ings sa 72 
yx 


The corresponding Eyring absolute rate equation 
for a bimolecular reaction in solution is!” 


moles\~! 
r ) , 8) 


liter 





RT YAYB 
kh =e—e—4Fexp/RTQAS*IR 
Nh x 


where AE.x, or Eact is obtained from the slope of 
log kvs 1/T and where we have added the ac- 
tivity coefficient term and imply a mole fraction 
equilibrium constant K,. 

Using the “cage model’’ of a liquid, Rabino- 
witch!*!° has developed an improved collision 
theory of bimolecular reaction rates in solution. 
His equation is 


Nw a*y; 
rerPaalee (: +—pe-# wt) : (9) 
No 4D; 


in which n is the coordination number (effective 
for reaction) of a reactant molecule, No is the 
number of liquid lattice points per liter, a is the 
‘lattice distance,” D; is an average diffusion 
coefficient of the reactants, v is a frequency 
characteristic of the solvent, y; is a factor ac- 
counting for the fact that a “‘jump”’ of molecule B 
may not destroy coordination with molecule 
A, and the other symbols have their usual 
significance. 

Equation (9) shares with earlier collision equa- 
tions based on gas kinetics the disadvantage of 
neglecting the entropy factor exp AS*/R arising 
from the internal degrees of freedom of the 
molecules concerned. It is nevertheless valuable 
because of the clear picture it gives of the 
significance of a “‘collision’’ and of the positional 
condition for reaction; namely, that two mole- 
cules must be in a state of coordination in order 
to react. Equation (9) also distinguishes clearly 





k=n 








17 Glasstone, Laidler, and Eyring, The Theory of Rate 
Processes (The McGraw-Hill Book Company, New York, 
1941), p. 199. 

18 Rabinowitch, Trans. Faraday Soc. 33, 1225 (1937); 
Rabinowitch and Wood, ibid, 32, 1381 (1936). 

#® Mark and Simha, Naturwiss. 25, 833 (1937). 


reactions with negligible energy of activation, 
and an observed dependence on viscosity from 
the commonly studied reactions with appreciable 
energy of activation and no observable depend- 
ence on viscosity. If Eact is small, the equation 
reduces to 


k=4nN,,D;/Noa?y: (10) 
and for larger values of Eact, to 


nNw 
k = ve-Fae RT, (11) 


No 





If No is equal to the total number of all mole- 
cules per liter and if m remains constant, the 
simple Eq. (11) predicts that k* will be inde- 
pendent of solvent. 

As it stands, Eq. (11) is evidently incompatible 
with Eqs. (5), (7), and (8). Arbitrary insertion 
into Eq. (11) of the missing entropy and activity 
coefficient factors, which must be found in any 
acceptable equation, gives 


YAYB 


Av 








yen Fact! RT eA St/R 


, (12) 


k=n 
No x 


Comparison of Eqs. (12) and (7) clarifies the 
significance of the mole fraction rate equation in 
terms of collisions. If the total number of lattice 
positions per unit volume, No, were a constant 
independent of the solvent, Eq. (12) would be 
essentially the usual volume molar rate equation. 
When however rate constants in solvents of 
differing molecular volume are compared in mole 
fraction units, the tacit assumption is made that 
the total number of molecules per liter present, 
the total number of positions in the liquid lattice, 
and the effective number of lattice positions for 
the reactant molecules are equal to each other. 
With these assumptions, Eq. (12) becomes essen- 
tially identical with Eq. (7). 


Temperature Coefficients of Ionic Reactions 


The energy and entropy of activation of an 
ionic reaction, calculated in the usual way at 
constant solvent composition, are functions of 
dD/dT for the particular solvent used.* 2?! It 


20 Svirbely and Warner, J. Am. Chem. Soc. 57, 1883 
(1935). 
21 Warner, Ann. N. Y. Acad. Sci. 39, 345 (1940). 
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has been pointed out that according to Eq. (1) 
the effect of electrostatic forces on the tempera- 
ture coefficient may be removed by extrapolating 
log k to infinite dielectric constant, 1/D=0. Thus 
values of the entropy, AS, and energy, AF}, of 
activation may be obtained, which are inde- 
pendent of the ionic properties of the reactants. 
Consideration of Eq. (1) suggests that the value 
of the entropy of activation so obtained is 
equivalent®*°*!_ to the so-called isodielectric 
entropy of activation, A.Stponst p obtained by arbi- 
traily varying the composition of the solvent in 
such a manner as to maintain the dielectric con- 
stant independent of temperature. 

Moreover, we note that by varying the solvent 
composition in such a way that the product DT 
is maintained constant, it is possible according to 
Eq. (1) to obtain a temperature coefficient and 
therefore an energy of activation, AF* const pr 
independent of the electrostatic terms and equal 
to AEt pu. 

In this case it is not necessary to compare 
results at zero ionic strength, but only at constant 
concentration of reactants, provided that the 
rate constant is separable into three terms as in 
Eq. (1) and that the ionic atmosphere term in- 
volves D and T only as the product; this is true 
not only of Eq. (1) but of all equations derived 
from the Debye-Hiickel postulates, such as those 
of Bjerrum” and of Gronwall, LaMer, and 
Sandved.” 

When the temperature coefficient is measured 
with D held constant, the entire electrostatic 
effect is arbitrarily thrown into the energy of 
activation, AE} onst p. Likewise when the product 
DT is held constant the electrostatic effect is 
arbitrarily absorbed in the entropy term, 
ASteonst DT. 

Equation (1) predicts that ASteonst p and 
AE} .onst pr Will be independent of dielectric con- 
stant and of solvent. Failure to find such an inde- 
pendence may be taken as indicative of the im- 
portance of solvent effects not accounted for by 
an equation derived from simple electrostatic 
postulates alone. In particular the consideration 
of AEt.onst pr May help to determine the relative 
importance of electrostatic and non-electrostatic 
effects even when a satisfactory extrapolation to 


22 Gronwall, LaMer, and Sandved, Physik. Zeits. 29, 558 
(1928). 
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zero ionic strength cannot be made. If ASteonst p 
and AEtionst pr really approximate closely the 
values for the reaction freed of electrostatic 
effects, they may be expected to be practically 
temperature independent. 


NEW EXPERIMENTS ON BROMACETATE- 
THIOSULFATE REACTION 


A. Materials 


Bromacetic acid was purified by fractional 
crystallization, sublimation under reduced pres- 
sure, and drying over sulfuric acid. Determi- 
nations of acidity and of bromine agreed with 
calculated values within 0.1 percent. 

C.P. methanol was distilled over sodium and 
silver nitrate. Concentration of solution was 
checked by density measurements. 

Glycine was purified by precipitation with 
alcohol from the aqueous solution, followed by 
three recrystallizations from water. A 1M solu- 
tion had a conductivity of 9 10-* mho and a pH 
of about 5.9. 

Atlas commercial d-mannitol was recrystallized 
twice from water. Melting point 165.9°. 

Dioxane and ?¢-butanol were purified as de- 
scribed by Tomlinson,® ¢-butanol receiving two 
redistillations. 

Ethylene glycol was stored for several days 
over calcium oxide and distilled under reduced 
pressure. 

Clear potato starch solution (about 0.1 percent) 
was used as an indicator in iodine titrations. The 
sensitivity of the end point was roughly 10-7 
mole per liter of iodine when D>70. 

Thallous chloride was purified and prepared as 
described by Hogge and Garrett." 

Analytical reagent potassium iodate, used as a 
primary standard, was recrystallized, ground and 
dried at 150°. C.P. sucrose and neutral potassium 
iodide were used. 


B. Methods 


The experimental methods are those previ- 
ously described by LaMer and Kamner* and 
Tomlinson® except as follows: 

Reactions in high dielectric solvents were 
stopped by dilution and cooling with a mixture 
of starch solution and chopped ice, followed by 
immediate titration from a weight burette. Re- 
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actions in solvents of dielectric constant lower 
than 40 were stopped by addition of a measured 
amount of iodine solution from a carefully cali- 
brated hypodermic syringe. The small excess of 
iodine was determined by measuring the trans- 
mittance of the solution at \=3650A with a 
Coleman spectrophotometer and comparing with 
a previously prepared curve, or by titration with 
dilute thiosulfate to 100 percent transmittance. 
The two methods gave identical results. 

When either dioxane or glycine was the added 
solvent, a sodium thiosulfate solution in this 
mixture, containing any added lanthanum chlo- 
ride, was mixed at zero time with a smaller 
quantity of a stronger sodium bromacetate solu- 
tion in water only, without added solvent. This 
precaution was taken to avoid peroxide formation 
in dioxane and hydrolysis of bromacetate ion in 
the presence of glycine. 

In glycine solutions it was important to add a 
little 0.1N acetic acid to the reaction mixture 
before titration. 

The method of Hogge and Garrett'® was em- 
ployed in following thallous chloride solubilities. 

The temperature of thermostats was de- 
termined to within 0.02° by means of a platinum 
resistance thermometer. 

Values of the dielectric constant of glycine 
solutions have been specially re-determined for us 
by Dr. Jeffries Wyman, Jr.** Other values have 
been taken from the literature.”* 


C. Results 


Data are presented in Table I for the brom- 
acetate-thiosulfate reaction in water, glycine- 
water, methanol-water, and sucrose-water at 25° 
and 45° and for /-butanol-water, glycol-water, 
and mannitol-water at 25°; in the lower dielectric 
range constants have been obtained in ¢-butanol- 
water, dioxane-water mixtures at 25° and 45° and 
in dioxane-water with added lanthanum chloride 
(Table II). In glycine solutions at 45° a small 
drift with time was observed, which was ap- 
parently caused by a glycine catalyzed hydrolysis 
of bromacetate ion. A short extrapolation to zero 
time, the uncertainity of which is less than 1 
percent, surmounted this difficulty. The drift can 


*3 \Vyman, in press. 
* Akerlof, J. Am. Chem. Soc. 54, 4125 (1932); Akerlof 
and Short, ibid. 58, 1241 (1936). 
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Solvent 


k log kx.o+1 1+ log ne, 


T=25°C a=6.4A 
0.1299 0.406 0.387 0.343 
.0815 -3220 355 327 
-0616 -2713 .302 .247 
.1446 3337 .212) 
.1238 308 313] 098 


% Sucrose 
MeOH 
7%, MeOH 
% MeOH 
MeOH 
MeOH 
% t-Butanol 
© t-Butanol 
Y t-Butanol 
% t-Butanol 
14.7% mannitol 
14.7% mannitol 
14.7% mannitol 
26.4% glycol 
2.309M glycine 
2.309M glycine 


.0914 .269 .213 
.0730 3174 356 .323 
-1460 4127 .323) 

.1333 .3934 .320} 

.1090 oa .322 

.1390 ‘ .405 

-1032 403} 

0945 : 404} 

.1660 46: .379 

-1666 1.021 .887 

extrapolated* 


T =45°C 
0.1701 
1722 
-2178 
-2639 
mean 
0.1991 
.1412 


a=6.8A 


2.406 
2.413 
2.677 
2.954 


Water 

Water 

Water 

Water 

Water 

0.363M glycine 
1.00M glycine 


3.101 
3.54 


.2447 
-1660 


4.25 
5.51 


1.00M glycine 

2.294M glycine 
glycine 

17.3% MeOH 


.1988 B: . 
31% sucrose a 


76- 
-1467 53 








* Compare Table III. 


be accounted for by postulation of a side reac- 
tion, second order in bromacetate ion and glycine, 
with a specific rate constant of 5 10~ liter per 
mole per minute. 

Recalculated data of LaMer and Kamner,' 
Straup and Cohn,‘ and Tomlinson® ° are given in 
Table III. The data of Kappana’ cannot be 
interpreted consistently and in particular his 
constants obtained in sucrose-water and ethanol- 
water mixtures cannot be extrapolated accu- 
rately to zero ionic strength, and are in disagree- 
ment with other work. They have, therefore, been 
omitted. 

In Table IV are listed a few new determina- 
tions of the solubility of thallous chloride, to- 
gether with some values from the literature. 


DISCUSSION 
A. The Bromacetate-Thiosulfate Reaction 


In Figs. 1 and 2 log ko for the bromacetate- 
thiosulfate reaction is plotted against the recip- 
rocal of the dielectric constant, with k calculated 
in 1 in the usual volume molar units, liters per 
mole per minute, and in 2 in the normalized mole 
fraction units of Eq. (3). The points obtained in 
glycine-water solutions lie on a straight line in 
both plots at 45° and 25°. When the solvents are 
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TABLE II. Summary of data obtained in low dielectric 
solvents. 








Dioxane-water 
25°C mixtures 


et M/l 
(min.~) log 555 


201.8 0.312 1.993 
135.9 275 1.858 
90.8 237 1.721 
58.6 .203 1.565 


t=45°C 
3.920 873 


3.907 878 
3.907 725 
3.895 363 


3.915 356 


No LaCls 
t =25°C 
4.034 0.702 
20.95 .895 
50.05 1.035 


1 =45°C 


2.98 0.301 
4.93 381 
3.96 357 
5.67 357 
5.09 339 


CLaCls =0.642 XCNa.S0; 


Concentration 
M/1-104 


3.963 
3.981 
3.984 
3.980 


log kN 





0.314 
373 
349 


20.1¢ 
23.3° 
25.3° 


3.958 
20.55 
20.17 
49.58 
49.58 


20.0¢ 
20.0° 
21.8° 
21.8° 
23.2° 








@ Dioxane water. + t-Butanol water. 


TABLE III. Summary of data recalculated and used for 
curves. T=25°C 








1 Hog kino 


0.398 
398 
.266 
.012 
508 
453 
506 


Solvent D 1 +log kx-0 


78.54 0.399 
78.2 401 


Observer 


K, F 





water 
1.47% sucrose 
32.4% 70.2 368 
67.9% glycerol 56.5 .286 
3.73M urea 88.4 555 
2.00M ‘“ 84.0 478 
400M “ 89.0 557 
5.00M “ 91.0 .584 520 
0.40M glycine 87.7 522 516 
062M “ 92.7 577 567 
100M “ 101.4 .669 655 
1.70M ethanol 74.0 .362 333 
341M “ 69.2 334 .276 
s1im 64.3 .308 .220 
2.60M glycine 138.1 .903 .869 
280M “ 142.6 917 .880 


SHNNNNNNNNMNRRRR 








K—Kamner, reference 3. 

S—Straup, reference 4. 

ete {LaMer and Fessenden, J. Am. Chem. Soc. 54, 2351 
(1932 

T—Tomlinson, reference 5. 


aqueous mixtures of dielectric constant lower 
than that of water, a linear plot may or may not 


be obtained: if methanol-water solutions are 
used, the points obey the linear relation of Eq. (1) 
within the experimental error ; if tertiary butanol 
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is the added solvent, there is a slight upward 
curvature. 

Comparison of Figs. 1 and 2 demonstrates that 
use of a mole fraction scale effectively eliminates 
most of the specific solvent differences evident in 
Fig. 1. The points obtained in water and in mix- 
tures of water with glycine, sucrose, mannitol, or 
methanol are collinear after normalization ; those 
obtained when the added solvent is ethanol, 
glycol, or glycerol show a slight deviation, while 
the points obtained in ¢-butanol-water mixtures 
lie appreciably above the others, but are much 
closer to the extrapolated high dielectric curve 
than when plotted on the volume-molar scale. 
The correction on the urea-points is too large, so 
that these lie below the common linear plot of 
Fig. 2. The values of the parameter, ra+rz, 
calculated from Figs. 1 and 2, (Table V) clearly 
illustrate the above conclusions. 

While isodielectric methanol-water and sucrose- 
water solutions are, when k” is used, equivalent 
solvents at 45° as at 25°, there remains a definite 
difference in the slopes of log kiLo vs. 1/D above 
and below the water point at the higher tempera- 
ture. The value of r4+72 obtained from the high 
dielectric curve does not change significantly be- 
tween 25° and 45°, while the value calculated 
from the sucrose-water points shows some tem- 
perature dependence. Such a temperature de- 
pendence suggests the existence of solvent effects 
not accounted for by Eq. (1) and not entirely 
removed by the use of k¥. The parameter ra+rz 
is a convenient catch-all for such effects. 

A linear relation between log ko and 1/D 
does not disprove the existence of non-elec- 
trostatic effects. 1/D is very nearly linear in the 
mole fraction of such solvents as methanol, 


TABLE IV. Solubility of thallous chloride. 








Solubility 
Solubility 55.3 
(moles/1) 


0.01617 
.01617 
01688 
.01530 
.01466 
.01373 
01405 
01893 
01213 


> moles/! 


0.01617 
01617 
.01869 
.01613 
.01494 
.01419 
01437 
.01923 
01217 


Solvent D 


H.O (reference 13) 78.54 
H,0 78.54 
16% d-mannitol 76.1 

8.3% glycol 76.1 

2.7% t-butanol 76.1 
5.5% methanol 76.1 
4.2% ethanol (reference 13) - 1 
0. 5M 
2 





glycine 9.99 
58 3 








* Kingerley and LaMer, J. Am. Chem. Soc. 63, 3256 (1941). 
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sucrose, glycol, ethanol, etc. (but not ¢-butanol) 
added to water, over a range from 0-40 percent or 
more. If there is a specific solvent effect, also 
linear in the mole fraction of added solvent, the 
plot of log ko vs. 1/D may still be straight. The 
values of log kp-. and r4+7z obtained will, how- 
ever, be in error and in general dependent on both 
solvent and temperature. The low and solvent- 
dependent values of r4+7rg obtained from the 
kinetics of carbinol formation of dye ions'*:*® are 
apparently at least partly caused by such a linear 
solvent effect. 

The data of Table I indicate that the Debye- 
Hiickel extrapolation to zero ionic strength can 
be made satisfactorily when D>60. It seems 
likely, therefore, that the deviations from line- 
arity in this range (Fig. 2) cannot be explained as 
the result of ion-association alone. 

Tables VI and VII contain values of the 
energies and entropies of activation and of the 
“frequency factor’ log B calculated from the 
Arrhenius equation and Eq. (8) using the data of 
Tables I, II and III and Figs. 2 and 3. A num- 
ber of conclusions may be drawn: 

In media of high dielectric constant the iso- 
dielectric B and AS? are practically independent 
of D, T, or solvent used, whereas the constant 
composition B and A.S* decrease as D decreases. 
This decrease corresponds to a lowered proba- 
bility of coordination of the like-charged reactant 
ions. 

AE? is temperature dependent when measured 
at constant composition, has a very slight de- 
pendence when measured in isodielectric media, 
and is practically independent of temperature 
when measured with the product DT held 
constant. 


TABLE V. Values of slope log kk=o vs. 1/D. 

















Added solvent Conc. scale Slope r4+rp (A) 
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Fic. 1. Volume molar rate constant k,—o (in min.) for 
bromacetate-thiosulfate reaction as function of tempera- 
ture and dielectric constant: 


ASt, calculated from the data for the lanthanum 
ion catalyzed reaction at a low constant ionic 
strength with dielectric constant. of 20 to 35, 
decreases in absolute value (becomes less nega- 
tive) as D decreases. This corresponds to a 
negative ionic strength effect consistent with 


/43 


400 


80 


60 


50 





pm. 


T 
—_ 


T 
Legend 


J-General linear curve 


Z-ethonol 





——e 


” 
#20 

> 

Ss 


N 


45° 


©-urea 
4-glycine 


IT-t-butanol 





k;, 


/+1og 





25 glycine, sucrose, Normalized —87.8 5.55 
methanol, mannite 
25 urea Normalized —69.1 7.04 
45 glycine Normalized —81.4 5.62 
45 sucrose, methanol Normalized —90.2 5.06 
0 sucrose Normalized —75.8 6.43 
25 mannitol vol. molar +18 —27 
25 sucrose vol. molar —20.6 +23.6 
25 t-butanol vol. molar —28.2 17.3 
25 methanol vol. molar —54.9 8.9 
25 glycine vol. molar —94.0 5.18 
45 sucrose vol. molar —35.3 12.9 
45 methanol vol. molar —59.5 7.7 
45 glycine vol. molar —91.3 5.34 
0 sucrose vol. molar —13.8 39 
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Fic. 2. Normalized mole fraction rate constant k)_o (in 


(194 en and LaMer, J. Am. Chem. Soc. 63, 3110 min.) for bromacetate-thiosulfate reaction as function of 


temperature and dielectric constant. 
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TABLE VI. Energies and entropies of activation in media 
of high dielectric constant. kY(min.~!) «=0. 











Solvent added —As+ 
Range to water D(25°) AEt logio B cal./deg. 
Constant Composition 

0-25° none 78.54 15,730 10.94 18.5 

0-25° 32.4% sucrose 70.2 15,530 10.64 19.9 
25-45° none 78.54 15,640 10.94 18.7 
25-45° 32.4% sucrose 70.2 14,850 10.19 22.1 
25-45° 17.3% methanol 70.5 14,980 10.26 21.8 
25-45° 1M glycine 101.4 15,660 11.12 17.8 
25-45° 2.3M glycine 131.3 15,850 11.47 16.2 

Constant Dielectric 

0-25° glycine 88.06 17 530 12.35 12.1 

0-25° sucrose 78.54 17,530 12.26 12.5 
25-45° glycine 78.54 17,720 12.38 12.1 
25-45° glycine 88.06 17,530 12.36 12.2 
25-45° sucrose-methanol 71.66 17,920 12.41 11.9 
25-45° t-butanol 71.66 17,270 11.97 14.0 
25-45° glycine inf. 15,960 12.23 12.8 

Constant (DT) 

0-25° glycine 80.67 16,180 11.30 16.9 

0-25° sucrose 71.9 16,030 11.05 18.0 
25-45° glycine 78.54 16,200 11.27 17.2 
25-45° sucrose 70.2 16,170 11.14 17.7 
25-45° glycine inf. 15,960 12.23 12.8 








previous conclusions®® that the reaction is be- 
tween the positive complex ion La+*+*-S,03> and 
the negative BrAc~. AF? const p is independent of 
dielectric constant, while AE? .onst pr decreases 
slightly as D decreases. This might be caused by a 
positive non-electrostatic effect and a corre- 
sponding dependence of parameters on tem- 
perature. 

When no lanthanum chloride is added a value 
of AEtwonst pr is obtained (Table VII) almost 
identical with that obtained for the catalyzed re- 
action. From this we conclude that the great in- 
crease in the reaction rate when lanthanum ion is 
present is purely an electrostatic effect. The value 
of log Beonst pr, On the other hand, is over two 
units higher for the catalyzed reaction. 

AE} wonst pt is about 1000 cal. lower in the low 
dielectric media than in the high dielectric media 
studied. This indicates that there may be an im- 
portant positive non-electrostatic effect in the 
dioxane-water and ¢-butanol-water mixtures and 
that the upward curvature of the plot of log k ws. 
1/D, found by Tomlinson® ® is not due entirely to 
the uncertainty of the extrapolation to zero ionic 
strength. 

If the value of ~Na/No in Eq. (12) is taken as 
1/55.5, and the entropy factor is taken as one, the 
value of the isodielectric frequency factor B 
calculated is only about 0.0-0.6 log unit 
larger than that experimently obtained, whereas 
Eyring’s estimates from the transition state 
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theory** for two reacting polyatomic molecules 
call for an entropy factor of 10-*-10-"". This ap- 
parent disagreement is typical of reactions in- 
volving ions and is presumably the result of large 
entropies of vaporization and large solvent and 
solvation effects. The agreement with predictions 
of the collision theory is probably quite fortui- 
tous; for many non-ionic reactions and for some 
ionic reactions (see, for example, reference 25) it 
is very poor. 


B. Effect of Choice of Concentration Scale upon 
Interpretation of Other Reaction Rates 


Search of the literature reveals some previous 
data, the interpretation of which might be 
affected by a change to the mole fraction concen- 
tration scale. Two interesting cases are discussed 
in Sections (C) and (D). Straup and Cohn‘ have 
studied the bromacetate-thiosulfate reaction in 
mixtures of water, ethanol, and either glycine or 
urea, iso-dielectric with water, with results con- 
sistent with those already discussed; the mole 
fraction constant is about the same in the iso- 
dielectric solvents, but appreciably low in the 
urea mixtures. 

For a reaction between an ion and a neutral 
molecule, the importance of the dielectric constant 
will depend upon the geometry of the critical 
complex. If the reactants can be considered as 
two particles approaching each other and forming 


TABLE VII. Energies and entropies of activation in media 
of low dielectric constant. k¥ in min.~!, t=25°—45°. 








Concen- + 
Solvent added tration of log B —AS+ 
to water reactants D AEt (min.)~! cal./deg. 





Constant Composition 
49.43% dioxane; LaCls 3.9 34.8 16,980 14.03 4.5 
61.09% dioxane; LaCls 3.9 25.0 15,720 13.38 7.8 
63.80% dioxane; LaCls 3.9 22.7 15,640 13.37 7.9 
66.50% t-butanol; 


No LaCls 20.5 23.3 15,960 11.29 17.1 
63.80% t-butanol; 

No LaCls 50.5 25.3 15,890 11.34 16.8 

Constant Dielectric 

dioxane; LaCls 3.9 20.2 13,960 12.21 12.8 
dioxane; LaCl; 3.9 22.3 14,000 12.19 12.9 
dioxane; LaCls 3.9 25.0 13,940 12.04 13.6 
dioxane; LaCls 3.9 30.1 13,940 11.92 14.2 
dioxane; No LaCls 3.9 20.0 13,810 9.66 24.5 
t-butanol; No LaCls 20.5 23.3 13,370 9.36 25.9 


Constant (DT) 


dioxane; LaCls 3.9 32.1 15,230 12.78 10.3 
dioxane; LaCl; 3.9 26.7 14,970 12.77 10.3 
dioxane; LaCls 3.9 23.8 14,990 12.86 9.9 
dioxane; LaCls 3.9 21.6 14,900 12.86 9,9 
t-butanol; No LaCls 20.5 23.3 14,480 10.18 22.2 
t-butanol; No LaCls 50.5 25.3 15,070 10.70 19.8 








26 Glasstone, Laidler, and Eyring, reference 17, p. 19. 
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Fic. 3. Effect of dielectric on La*** catalyzed reaction. 


a double-sphere critical complex such as that of 
Eq. (1), with the charge continuing to be cen- 
tered on the reacting ion, little effect will be ex- 
pected. This is an extreme case. The opposite 
extreme, in which the critical complex can be 
represented as a single sphere with the charge 
symmetrically distributed about its center leads 
to the equation?’ 


Zare*Nn Zare?Nw 
2DRTr, 2DRTrx 





In kyo =In Rp» + (13) 


Since r, is greater than 74 Laidler and Eyring con- 
clude that the rate of a reaction between an ion 
and a neutral molecule should increase with 
decreasing dielectric constant. 

As the actual case probably lies somewhere 
between these extremes, Eq. (13) should predict a 
self-potential effect of the right sign but too large. 
The experimental examples suggested?” *8 are not 
entirely convincing. The supposed dielectric con- 
stant effect on the reactions between water and a 
sulfonium cation” and between alkyl halides and 
halide ions*® is so large that one suspects other 
effects may be very important. On the other 
hand, the reaction between hydroxy] ion and iso- 
propyl bromide, studied*! in 60 percent and 80 
percent ethanol solutions is significantly faster in 

27 Laidler and Eyring, Ann. N. Y. Acad. Sci. 39, 303 


(1940). 
28 Glasstone, Laidler, and Eyring, reference 17, p. 439 


q. 
sag and Ingold, J. Chem. Soc. (London), 1571 
30 Bergmann, Polanyi, and Szabo, Zeits. f. physik. Chemie 
B20, 161 (1933). 
3t Hughes, Ingold, and Shapiro, J. Chem. Soc. (London), 
225 (1936). 


se 


the latter solvent only when volume molar con- 
centration units are used. On a mole fraction 
basis the “‘effect’’ will evidently become of the 
wrong sign at a temperature slightly below those 
used for the measurements (Table VIII). Proba- 
bly no definite conclusions as to the importance 
of dielectric constant effects should be drawn 
until it is possible to separate the non-electro- 
static effects. 

There are, however, at least two similar cases 
in which this can be accomplished in part. The 
increase of the rate of reaction between thio- 
sulfate ion and methyl bromacetate®*® when 
t-butanol is added to the water-solvent, is more 
than removed by conversion to a mole frac- 
tion scale. Comparison with the observed sol- 
vent effects on the corresponding reaction with 
bromacetate ion suggests that in this case the 
self-potential effect is unimportant. The reaction 
between alkyl halides and thiosulfate ion has 
been studied‘ in water, in mixtures of water with 
ethanol, glycine, and urea, and in mixtures of 
water, ethanol, and either glycine or urea iso- 
dielectric with water (Table IX). The results 
indicate a small self-potential effect of the sign 
predicted by Eq. (12) when glycine or urea is 
added, although this is only evident in the latter 
case when the mole fraction scale is used. When 
ethanol is the added solvent a similar effect 
appears to be present, but is far less important 
than another solvent effect. 


C. The Quenching of Fluorescence 


The quantity (Io/I—1)/Cg=ke is obtained® 
from studies of the quenching of fluorescence in 
solution. J) is the intensity of fluorescent light 
emitted by a solution before addition of a 


TaBLE VIII. Effect of solvent upea reaction between 
isopropyl bromide and hydrexyl ion, from data of 
Hughes, Ingold, and Shapiro, reference 31. 











kN (60% 
4 Medium D k ethanol) 
80 80% Ethanol 23.20 4.61 3.52 
80 60% zi 31.82 3.11 3.11 
50 80% i 28.10 0.241 0.184 
50 60% as 51.48 181 181 
35 80% = 30.8 .0477* .0345* 
35 60% = 41.1 0368" 0368* 








* Values obtained by extrapolation. 
® Wawilow and Frank, Zeits. f. Physik 69, 100 (1931). 
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TABLE IX. Effect of solvent upon reaction between 
methyl iodide and thiosulfate ion, from data of Straup 
yo, ae reference 4. CH;I1+S20;-—-CH;S.0;-+I-. 
t= 25°C. 








Medium D logk—logky.9 — log RN —log ky. 





0.241 0.183 


.204 141 
054 060 
395 307 


.241 145 
74 75 


353 21 
022 042 


69.5 


78.5 
87.7 
64.5 


78.5 
92.7 


78.5 
91.0 


20% ethanol 
0.4M glycine in 
20% ethanol 
0.4M glycine 
30% ethanol 
0.62.M glycine in 
30% ethanol 
0.62M glycine 
5M urea in 30% 
ethanol 
5M urea 


i> +1 


++ 








quenching material and J is the intensity at a 
concentration Cg of quencher. A constant source 
of illumination is used, and the molar concen- 
tration of fluorescent material is in general much 
smaller than that of quencher. To show that ka 
is proportional to the second-order rate constant 
for the quenching reaction, we need make no 
assumptions concerning the mechanism*®—** be- 
yond that the process of activation of a molecule 
F by a light source, the emission of fluorescent 
light, and the deactivation of an activated mole- 
cule by a quenching molecule may be represented 


by equations of the following type: 

ky 
F+hv-F* zero order for constant 
Cr and constant light 
source, 


ke 
Ft F+hyp 
k 
F t4+0—F +0Q+Energy second order in Ft and Q. 


first order in F?, 


Since J is independent of time a steady state 
exists. 


dCrt/dt=0=ki—koCr*—k3CriCg. (14) 


We choose the units of k; and k2 so that J and Jo, 
in Einsteins per second, are given by 
Io=hi, 
IT=k2Crt. 
Substituting Cr? from Eq. (14), we have 


Io 1) /c ( ky 1) /c ks k 
z oe NeeCet — 


% Weiss, Trans. Faraday Soc. 35, 48 (1939). 

* Bowen, Trans. Faraday Soc. 35, 15 (1939). 
( % Rollefson and Stoughton, J. Am. Chem. Soc. 63, 1517 
1941). 


(15) 
(16) 


(17) 
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To obtain k; the true rate constant of the 
quenching reaction from kg, it is necessary to 
multiply by ke. This is equivalent to dividing by 
the mean life of fluorescence, a quantity which 
may be determined for any particular fluorescent 
substance by direct optical measurement,*® or 
indirectly by studying the effect of viscosity on 
the polarization of the fluorescent emission.*” ** 
Since quenching has an activation energy of 
about 4 kcal., i.e., similar to that of diffusion, a 
linear relation between kg and diffusion (fluidity) 
has been postulated by many authors,*® #8 but 
experimental curves show deviations from line- 
arity (see Fig. 4.1) which have not received 
satisfactory explanation. 

When the reactants are ions, the almost obvi- 
ous dependence upon ionic strength has been 
observed.‘ *! There is also a self-potential effect. 
In this connection it should be pointed out that 
any attempt*® to base a decision of the mecha- 
nism between the concept of “‘an ordinary bi- 
molecular reaction between the quenching ion 
and the activated dye molecule’ and that of 
“de-activation of the photo-activated molecule 
by collision with the quenching ion in which some 
energy is transferred to the quencher’”’ by testing 
the agreement with the Brénsted activity coeffi- 
cient factor is doomed to failure. The number of 
effective collisions is not determined by the ac- 
tivities of the colliding molecules alone. Brénsted’s 
equation can be justified equally well from an 
equilibrium viewpoint as from a purely collisional 
or diffusional viewpoint. (Reference 3, p. 2668.) 

Since most investigators of fluorescence quench- 
ing have been unaware of the importance of these 
electrostatic complications, there are few data in 
the literature that can be interpreted unambigu- 
ously. It is evident, however, that the plot of kg 
vs. fluidity (1/n) is never linear, whether the 
viscosity (n) is varied by changing the medium, 
or by varying the temperature. Figures 4 to 7 
demonstrate that when the medium is varied at 
constant 7, use of the normalized mole fraction 


36 Gaviola, Zeits. f. Physik 42, 853, 862 (1927); Szyma- 
nowski, ibid. 95, 440 (1935). 
37 F, Perrin, Thesis, Paris, France (1929). 
38 Wawilow, Acta. Physicochim. U.S.S.R. 7, 49 (1937). 
89 SveSnikov, Acta Physicochim. U.S.S.R. 7, 755 (1937). 
4° Stoughton and Rollefson, J. Am. Chem. Soc. 61, 2634 
(1939). 
( on and Rollefson, J. Am. Chem. Soc. 62, 2264 
1940). 
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constant kg” reduces the curvature considerably.* 

The quenching of fluorescein® (Fig. 7) and of 
rhodamine-B® (Fig. 6) by potassium iodide are 
ionic reactions. The most interesting result in the 
latter case is that obtained when pure methanol 
is the solvent. In this case kg is found to be much 
greater than in water, although iodide ion is 
known to diffuse more slowly in methanol. After 
normalization the ratio of kg in water and in 
methanol is 1.40, while the ratio of the mobilities 
of iodide ion is 1.25 and of the diffusion constants 
of potassium iodide is 1.26. 

The quenching of acridone by potassium 
iodide“ and the quenching of rhodulin red by 
aniline® involve a neutral molecule and should be 
much less dependent upon dielectric constant. In 
both cases the plot of ke” vs. 1/n is linear. The 
case of rhodulin red is particularly interesting, 
since the quencher is the comparatively large and 
neutral aniline molecule, which might be ex- 
pected to obey Stokes’ law, and since pure 
alcohols were used as solvents. 

It can hardly be a coincidence that the percent 
deviation from a straight line through the origin 
of points, obtained in a series of alcohols of 
molecular weights ranging from 32 to 85, should 
be closely proportional to the molecular weights. 
It is to be noted that the anomaly in the order of 
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Fic. 4, The quenching of fluorescence of rhodulin red by 
aniline. Data of SveSnikov. 


* The symbol kg or kg’, used in the figures, denotes the 
first-order constant (Jo/J—1) obtained with Cg held 
constant. 
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Fic. 5. The quenching of fluorescence of acridone by potas- 
sium iodide. Data of Stoughton and Rollefson. 


the points obtained in isobutyl and isoamyl 
alcohols is removed by normalization to the mole 
fraction scale. 

The significance of this linear plot is not 
entirely clear, and cannot be explained, at least 
until a thorough experimental analysis of the 
problem, involving a determination of diffusion 
coefficients, has been made. It must be empha- 
sized that a mere change of scale cannot reduce 
the curvature of the plot of k vs. 1/n when the 
latter is varied by changing the temperature and 
holding the solvent composition constant; and 
in the case of the quenching of rhodulin red by 
aniline, certain points obtained in glycerol- 
ethanol mixtures fall on a smooth curve with the 
points obtained in pure solvents before normaliza- 
tion to a mole fraction scale, but apparently do 
not fall on the linear plot after normalization. 

Nevertheless the quenching of fluorescence is 
an example of a type of reaction in the study of 
which the concentration scale is very important 
and for which we find use of the mole fraction 
scale gives a simple result. Quenching involves a 
diffusion process as a rate determining step. Ex- 
periments are in progress to determine, if possi- 
ble, whether the cage model Eq. (10) or the 
Smoluchowski type equations fit the data best. 
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D. Other Solvent Effects. The Solubility of 
Thallous Chloride 


The rate of isomerization of ammonium 
cyanate!>” and of the carbinol formation of 
bromphenol blue" is slower in glycol-water 
mixtures than in the isodielectric methanol-water 
solutions. In such a case it is evident that use of 
kN will make the discrepancy worse. The order, 
by effect on the rate, of the added solvents, 
methanol, ethanol, glycol, is the order of their 
increasing molecular weights and viscosities ;!5 
but whether we reason from an equilibrium rate 
theory or from Eq. (9) we must conclude that 
viscosity cannot directly affect the rate, since 
large energies of activation are involved. 

If the equilibrium theory of reaction rate is 
correct, solvent effects should not be peculiar to 
ionic kinetics, but should be found for ionic 
equilibria in general. We have, therefore, meas- 
ured the solubility of thallous chloride in several 
isodielectric solvents (Table IV). While no such 
extended analysis as that given by Hogge and 
Garrett for ethanol-water mixtures has been 
attempted, it is possible to draw certain con- 
clusions: 

The solubility cannot be predicted unambigu- 
ously from the dielectric constant. It depends 
upon the solvent added to water and increases in 
the order of added solvents methanol, ethanol, 
t-butanol, glycol, mannitol. This is consistent 
with a decrease in reaction rates, which involve 
the factor fafe/fx. A general decrease in the 
activity coefficients, corresponding to a _ pro- 
portional increase in the solubility of reactants 
and critical complex, would decrease the reaction 
rate. The solubility of thallous chloride in water 
is actually increased by adding mannitol. Dielec- 
tric constant effects are here of less importance 
than another effect, apparently related to the 
size of the mannitol molecule. The addition of 
glycine increases the solubility, but less than 
would be predicted from the ethanol-water data, 


CONCLUSIONS 


One conclusion may be drawn; namely, that 
use of the mole fraction scale instead of the 
volume molar scale leads to distinctly simpler 


4 Svirbely and Schramm, J. Am. Chem. Soc. 60, 330 
(1938); Lander and Svirbely, ibid. 60, 1613 (1938). 
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and more consistent results in various solvents in 
the specific case of the bromacetate-thiosulfate 
reaction. 

The more general conclusion, that the mole 
fraction scale should always be used in com- 
paring rate constants in different solvents, is in a 
sense an academic one requiring no experimental 
verification. As Scatchard clearly states! the only 
object is to choose that scale which will allow the 
specific solvent effects to be most readily ac- 
counted for by activity coefficients determined 
by the physical properties of solute and solvent. 

The basis of any mole fraction ideal law of 
solution is the assumption that the number of 
available positions in a liquid is proportional to 
the total number of molecules present and that 
all positions are equally available for all mole- 
cules, regardless of size or shape. Then any size- 
shape effects must be compensated for by activity 
coefficients. If there are other deviations from 
ideality, i.e., if the forces between solvent- 
solvent, solvent-solute, and solute-solute pairs of 
molecules are appreciably different, these devia- 
tions must also be compensated by an activity 
coefficient. 
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Fic. 6. The quenching of fluorescence of rhodamine B by 
potassium iodide. Data from Frank and Wawilow. 
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Fic. 7. The quenching of fluorescence of fluorescein by 
potassium iodide. Data from SveSnikov. 


There are no equations available for precise 
calculation of non-electrostatic activity coeffi- 
cients. The approximate equation valid, for 
“regular”’ solutions?’ 


a 


En , Exe a 
()-Ga fo 
V1 V2 


where the subscripts 1 and 2 refer to solute and 
solvent, respectively, V = mole fraction, v= molar 
volume, and E’s are potential energy terms that 
may be replaced by the energies of vaporization 
of the pure liquids, predicts a relationship be- 
tween solvent effects and the ratio E22/v2. This 
however does not vary significantly in the aque- 
ous mixtures discussed in this paper. 

Except in certain special cases the effect of 
differences in size of solvent and solute molecules 
cannot be neglected. In the case of approximately 
rod-shaped molecules it has been shown theo- 
retically and experimentally*“ to be zero: Thus, 
n-hexane and n-hexadecane obey the mole frac- 
tion ideal laws throughout the entire range of 
concentration. For other shapes, it has been sug- 


48 Hildebrand, J. Phys. Chem. 43, 109 (1939) ; Hildebrand 
and Sweney, ibid. 43, 297 (1939). 

“* Hildebrand, J. Am. Chem. Soc. 59, 794 (1937). 
ass and Rushbrooke, Trans. Faraday Soc. 33, 1272 
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gested*®4® that there may be large deviations 
from the mole fraction laws and that these 
deviations are an increasing function of the ratio 
of the molar volumes of solute and solvent. For 
spherical molecules A and B with v4=1/2vz, the 
equations of Chang lead to a value of about 0.8 
for the activity coefficient of the very dilute 
solute, A.*7 

The simplicity of the results found for the 
bromacetate-thiosulfate reaction must be due toa 
considerable cancellation of the non-electrostatic 
activity coefficients of the reactants with those of 
the critical complex. No such cancellation is 
likely in the case of solubility or for a reaction 
such as the carbinol formation of brom-phenol 
blue“ in which the critical complex closely 
resembles one of the reactants. 


SUMMARY 


New data are presented to emphasize the 
difficulties of comparing reaction rate constants 
obtained in different solvents. It is shown that 
use of a mole fraction concentration scale instead 
of the usual volume molar scale gives a simpler 
result for the kinetics of the bromacetate- 
thiosulfate reaction. 

Temperature coefficients have been determined 
for this reaction and are discussed as criteria of 
electrostatic and non-electrostatic effects. 

It issuggested that determination of the energy 
and entropy of activation with the product DT 
and the concentration of reactants held constant 
helps in interpreting the kinetics of ionic 
reactions. 

Data on the solubility of thallous chloride in 
solvents of different molecular volume are pre- 
sented to show that there are other specific 
solvent effects as important as the electrostatic 
effects. 

Examples of the effect of choice of concentra- 
tion scale upon the interpretation of data are 
taken from the literature. The quenching of 
fluorescence in solution as determined by the 
viscosity is discussed. 


46 Chang, Proc. Camb. Phil. Soc. 35, 265 (1939). 
47 Fowler and Guggenheim, Statistical Thermodynamics, 
(Cambridge University Press, 1939) p. 366. 
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An equation has been derived for the reaction rate constant in the case of ion-dipole reactions. 
Section II contains the extension of Debye-Hiickel’s theory of ionic atmospheres to the case of 
dipoles surrounded by ions, coupling it with Onsager’s theory of electric moments of molecules in 
liquids. Section III applies the potential constructed to the derivation of the rate constant by the 
kinetic procedure of Christiansen and Scatchard. The result is substantiated by existing data 
both from concentration and dielectric constant dependence of the rate of reaction between ions 


and dipoles (Section IV). 





I. INTRODUCTION 


INCE the formulation of the law of mass 
action and the promulgation of the theory of 
temperature coefficients, the two outstanding 
contributions to the theory of the kinetics of 
reactions in solution have been the application 
by Brénsted! of the theory of Debye and Hiickel 
to reaction rates, and Scatchard’s? development 
of the effect upon reaction velocities of changing 
solvent insofar as this effect is due to changing 
dielectric constant. These two contributions were 
applicable to the case of reactions between ions 
and led to the general kinetic equation 


ZaZpe? ZAZpe- K 


In k=In ko ’ 
DkTro DkT 1+ax 





(1) 


which has become so familiar to those who read 
kinetic literature. 

The authors, studying the reactions between 
charged ions and dipolar molecules, have become 
convinced that a similar equation expressing the 
concentration and dielectric effects upon ion- 
dipole reactions in solution should be available 
to kineticists. This paper has as its object, there- 
fore, the derivation of such an equation together 
with its experimental verification insofar as 
existing data permit. 


Il, THE IONIC ATMOSPHERE AROUND A 
DIPOLAR MOLECULE 


In determining the ionic atmosphere around a 
dipolar molecule we shall disregard the inter- 
1J. N. Bronsted, Zeits. f. physik. Chemie 102, 169 (1922); 


115, 337 (1925). 
2G, Scatchard, Chem. Rev. 10, 229 (1932). 


action between dipoles. Under this restriction the 
differential equation for the potential Yo in the 
neighborhood of any one dipolar molecule will be 
the same as for the potential around an ion. It is 
only with regard to the boundary conditions at 
the surface of the molecule that the two problems 
differ. 

We shall assume that there are s different 
species of ions. Let 2; be the valency of the species 
z and let n; be their number per cm*. Then the 
differential equation for Yo becomes, with the 
usual restriction to small values of Wo, 


Vpo= Ko, (2) 


where « is the well-known parameter of the 
theory of Debye-Hiickel* and is of the dimension 
of a reciprocal length: 


x2 = (4re?/DkT) > n32;?. (3) 
i=1 

Here ¢ signifies the elementary charge and D the 

dielectric constant of the solvent. 

As the Eq. (2) can be ‘‘separated”’ in spherical 
coordinates it is easy to obtain a general solution 
which can be adapted to any distribution of 
point charges within the molecule. However, we 
shall be satisfied to treat the dipole as a point 
singularity. Then a suitable particular solution 
of (2) can be obtained in a very direct way by 
“‘polarizing’’ the solution of Debye and Hiickel. 
This is, for the case of an ion, 


Yo=Ace""/r. (4) 


3 P. Debye and E. Hiickel, Physik. Zeits. 24, 185 (1923). 
See also Hans Falkenhagen, Elektrolytes, English transla- 
tion by R. P. Bell (Clarendon Press, Oxford, 1934), 
Chapter VII. 


598 





mn fn fy -—-= © 


_- ~ DD &¢ 


KINETICS OF 


Let the position of the dipole coincide with the 
origin and its direction with the positive z axis. 
Let, furthermore, the shape of the dipole-mole- 
cule be idealized by a sphere of radius a. Then we 
obtain an adequate solution for r=a, i.e., a 
solution that has the same angular dependence 
as a dipole-potential, by differentiating (4) par- 
tially with regard to z. Doing this and sub- 
stituting in the result for z the expression 
rcos? (the angle 3 being counted from the 
positive z axis) we obtain 


Yo=A(e-"/r’?)(1+kxr) cosd, r=a. (5) 


As for the interior of the dipole molecule we 
might follow the procedure of Debye and Hiickel 
by determining y; there from the differential 
equation 


Vy¥i=0, rSa (6) 


and retaining the constant value D also for the 
interior. The condition to be imposed on y in 
this domain is that it should be continuous except 
for a singularity at the origin corresponding to a 
dipole of the given strength. 

If this is done a reasonable formula for the 
ionic atmosphere is obtained. However, a com- 
parison with experimental kinetic data (see 
Section IV) shows that, in order to obtain 
agreement, the dipoles have to be chosen of the 
order 10X107'8, i.e., five to ten times too large. 
The reason for this discrepancy is evidently that 
in the indicated procedure the interaction be- 
tween the dipole and the dielectric solvent is not 
properly accounted for. 

From Onsager’s theory of electric moments of 
molecules in liquids‘ it follows that any one 
immersed dipole acts upon distant charges with 
an ‘‘external moment”’ different from its moment 
in vacuum. The intervening medium has the 
effect of ‘enhancing both the permanent and 
the induced dipole moments.”’ We are, therefore, 
going to combine Debye-Hiickel’s theory of 
ionic atmospheres with Onsager’s model of a 
dipole immersed in a dielectric liquid. 

In accordance with Onsager we retain the 
spherical shape of the molecule and characterize 
it by its permanent moment yo (in vacuum) and 
its polarizability a. This is to be related to an 


*L. Onsager, J. Am. Chem. Soc. 58, 1486 (1936). 
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“internal refractive index n’’ by 
a= ((n?—1)/(n?+2))a’. (7) 


The problem to be solved is the following: a 
solution of (6) having a dipole singularity at r=0 
has to be joined to our solution (5) in the proper 
way (see reference 4, p. 1488) at r=a. This in- 
volves a determination of the external charac- 
teristic moment yu* such as it becomes under the 
action of the polarization created by itself. 

First we solve the problem which presents 
itself if a rigid dipole of moment m is introduced 
into a cavity of radius a. For the exterior the 
solution given by (5) is to hold; for the interior 
the solution 


Yi:=(m/r+Br) cosd, rsa (8) 


meets the requirements. From the usual boundary 
conditions for the electrostatic field the constants 
A and B are determined as ; 


A =3me"*/(D(2+2xa+«’a?)+(1+«a)) (9a) 
and 


B= —(m/a*)(D(2+2xa+«a*) —2(1+xa))/ 


(D(2+2xa+«*a*)+(1+«a)). (9b) 


Now m remains to be determined as a function 
of uo and n?. This determination has to be based 
on the “condition for internal equilibrium” 
(Onsager, reference 4, p. 1488): 


m=potaF., (10) 


where the force F, is the force acting on the 
dipole owing to the polarization created by itself. 
This force is found to be 


fe) 
F,= ——(¥;—m cos 3/r*?) = —B. 
02 


(11) 


Making use of (7) and (9b) we obtain 


m= po(n?+2)(D(2+2Kxa+ xa") + (1+ x«a))/ 


3(D(2+2xa+«a?)+n?(1+«a)). (12) 


With this value introduced into (9a) the poten- 
tial of the ionic atmosphere around a dipolar 
molecule of moment yo can finally be written in 


the form 


Yo= (u*e-*"/ Dr’) (1+ xr) cos 8, (13) 
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where u*, the ‘‘external moment,” is given by 


u* = wo(n?+2)De*/(D(2+2xa+x*a") 
+n*(1+«a)). (14) 


In the absence of an ionic atmosphere, i.e., 
for x0, this solution and the corresponding one 
for ¥; reduce to Onsager’s solution, as they 
should. For n?—1 Eqs. (13) and (14) represent 
the simpler case where the molecule has a 
permanent dipole but no polarizability. 


Ill. THE KINETIC EQUATION 


We shall now apply the potential constructed 
in Section II to calculate the rate of formation 
of an intermediate complex X formed out of a 
dipole molecule A of moment yo and an ion B of 
charge ezg. In doing so we are going to follow 
the procedure of Christiansen and Scatchard. 
However, it must be taken into account that 
the potential does not depend on r only but 
also on #. 

Let C.,, Ca, and Cg be the concentrations of 
X, A, and B, respectively, and let us introduce 
polar coordinates 7, 3, g about the center of a 
dipole A. The probability of finding an ion B in 
a specified element of volume characterized by 
the limits 7 and r+dr, 3 and 3+d#, g and g+de¢, 
will then, by Boltzmann’s theorem, be propor- 
tional to 


CaCp exp (— oeze/kT)r’ sin ddrdddg. (15) 


In order to obtain the rate of formation of X, 
and thereby the velocity of reaction, some 
criterion for the formation of X must be estab- 
lished. This can be done in two different ways. 
Either it is assumed that there are ‘‘sensitive 
zones’ on the exterior of the molecule. Whenever 
one of these areas is touched by an ion B the 
formation of X ensues. Or it is assumed that 
there is a critical distance for each direction. 

Let us follow the first alternative. This amounts 
to assuming that 7, 3, g have to assume specified 
values between comparatively narrow limits: 
ro and ro+Aro, Fo and Io+Ado, go and gotAgo. 
We shall assume rotational symmetry of the 
molecule about the z axis. Then ¢ need not be 
specified and (15) can be integrated with regard 


5J. A. Christiansen, Zeits. f. Physik. Chemie 113, 35 
(1924); G. Scatchard, J. Am. Chem. Soc. 52, 52 (1930); 
Chem. Rev. 10, 229 (1932). 
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to y. Consequently the expression for the concen- 
tration of the complex becomes 


C.= Qrk'’CaCare? sin VoAr AV 
Xexp (—ezey~o/kT), (16) 


where k’ is a constant.® 
Setting 
K= 2rk' ry? sin ToArpAVo (17) 
we obtain 


In (C,/CaCpg) =|In K—ezppo(ro, 3o)/RT. (18) 
In the limiting case x0 this reduces to 


In (C°/Ca°Cp°) =In K 
as €Zpuo™ COs 3o/Dk Tr, ( 19) 


where yo* is Onsager’s value 
bo* = wo(n?+2)D/(2D+n’). (20) 


From (18) and (19) we finally obtain for the 
speed of formation 


In R=In k,~0 
€Z~ COS Uo 


+ 
DkTr,? 





(uo* —p*e—*"°/(1+«ro)). (21) 


This formula for ion-dipole reactions corresponds 
to formula (1) in the case of ion-ion reactions. 

In the second alternative Eq. (15) has to be 
integrated to give the resultant probability of 
formation. In the most general case the critical 
distance will be a function of the angles # and ¢ 
and this functional relation constitutes the ‘‘shape 
of the molecular dipole.” It is, however, required 
that 7 is, in all directions, larger than a as other- 
wise Yo is not the valid potential throughout. 
The further procedure, after (15) has been in- 
tegrated follows the same lines as above. 

It is easy to perform the calculations if 7 is 
assumed constant, i.e., between 79 and ro+Aro, 
and if # is unrestricted. This means that the ion 
has to approach, in all directions, the same 
critical distance ro. Such a behavior would be 
rather astonishing for a molecular dipole. As a 
matter of fact, this assumption leads to the 
wrong sign for the dependence of k on concentra- 
tion and dielectric constant. 


6 It should be mentioned that, according to the kinetic 
theory, k’ is still proportional to 7! but this is of no im- 
portance for what follows. 
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A more adequate assumption would be that 
the critical surface is a prolonged spheroid with 
the dipole in one of the foci directed along the 
long axis. We have carried through the calcula- 
tions for this case restricting ourselves to weak 
potentials (as in the derivation of Debye’s Eq. 
(2)). The result is of a similar nature as (21) 
and coincides, for small values of «, with (21). 
The only difference which remains is that cos # 
is replaced by e, the eccentricity of the spheroid. 
This shows that the result may be chosen with 
the required sign according to whether the posi- 
tive or the negative side of the dipole is nearer 
the surface. 

We do not want to stress these calculations 
which involve a spheroidal form for the following 
reason. A spheroid would partake of the heat 
motion with rotations and these would have such 
a high frequency that the ionic atmosphere 
would not be able to follow them owing to its 
considerable time of relaxation. 

However, these’ calculations show that it is 
legitimate to choose 3» and fo as critical average 
values which are characteristic of the molecule. 
In this sense we are going to understand them 
and, consequently, use formula (21). 

It will be seen in Section IV that cos > has 
to be chosen as +1 if we want to obtain agree- 
ment with empirical data for positive ions, i.e., 
an increase of In k with an increase of 1/D. 
This means, in the spheroidal model, that the 
negative charge of the dipole has to be nearer 
the surface if the ion is positive. Furthermore, 
we assume that the critical distance ro coincides 
with the radius a used in Section IT. 

Introducing the dimensionless variable 


Z=ka=kPo (22) 
and the dimensionless quantity 
W= (In R—In Ryo) (2DRT ro?) /ezpuo* cos 8p = (23) 
we can write our result in the final form 
W=2?/(1+2+2?/2+(n?/2D)(1+2)). (24) 


This formula brings out the dependence of k on 
concentration in a general way. Whatever the 
valency of the B ion, and whatever the nature 
of the other (admixed) ions contributing to the 
ionic atmosphere, the relationship between W 
and z ought to be the same. 


As for the dependence of the reaction rate on 
D, the dielectric constant of the solvent, an 
additional remark should be made. Starting from 
the idea that k,o, as used above, may still 
depend on the value of D, several authors’ have 
split off from k,.o the part k,.9 measuring thus 

D=0 

the D dependence of k from a point of reference 
defined by the double transition x=0 and D- 
(see Eq. (1)). This procedure has proved to be 
successful in the case of ion-ion reactions. How- 
ever, in our case it would lead to the result that 
the change of k with 1/D should be in the 
opposite sense to the change with x. This is 
contradicted by all available empirical data. 

It should be taken into account that, in our 
case, the situation is very different from what it 
is in ion-ion reactions. In the latter case the 
electric attractions tend towards zero as D in- 
creases indefinitely. However, for a dipolar mole- 
cule, the ‘‘external moment’’ increases with in- 


crease of D. Therefore, though the lim y* exists, 
D=2 


the point D=~« is certainly not an adequate 
point of reference since it accentuates rather than 
eliminates the effect of the dielectric constant. 

In the case of ion-dipole reactions, owing to 
the argument just presented, the point «=0 can 
be interpreted as the ideal limiting case in which 
all electrostatic actions have disappeared. For 
this reason it seems to us that this point repre- 
sents the adequate point of reference. Anyway 
we are going to introduce this as a special assump- 
tion. Then our formula (21) will express the 
complete dependence of k on D as on concen- 
tration. 

Again it is useful to introduce dimensionless 
variables. If we set 


=)?/D (25) 


will still be of the dimension cm=', but free 
from D. Setting then 


c= Aa=Aro (26) 
and 


W’ = (In R—In ky~0)(2RkT)/ezpuo*d? cos Bo (27) 
we have the general relation 

W! =(1/D2)(1/(1+¢/D!+¢2/2D)). (28) 

7V. K. LaMer, J. Frank. Inst. 225, 709 (1938); E. S. 


Amis and V. K. LaMer, J. Am. Chem. Soc. 61, 905 (1939); 
J. C. Warner, Ann. N. Y. Acad. Sci. 39, 345 (1940). 
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It shows that, in the limit A-0, W’ will 
change as 1/D*. For larger concentrations the 
increase will become less and depends on the 
parameter ¢. It should be mentioned that 2? as 
well as ¢° are, by their definitions (3) and (25), 
proportional to the ionic strength. 


IV. APPLICATIONS 


Equation (21) may be tested as to the effect 
of ionic strength upon reaction rate by putting 
it into the non-dimensional form given by Eq. 
(24) and plotting a theoretical curve W vs. 2°. 
As a good approximation the term (n?/2D)(1+32) 
may be omitted in making the plot. In Fig. 1 
the solid line is the theoretical curve for n?=0 
and the broken line that for ?=9. It can be 
seen that for these extreme values of n? the plots 
do not vary a great amount. 

Experimental values of (In k—In k,~o) plotted 
against ionic strength will confirm the theoretical 
curve if they can be made to fit it by change 
of scale, for 2 is proportional to the ionic 
strength. From the change of scale of the 
abscissae a reasonable value of ro ought to 
result. Further, from two given values of log k, 
corresponding to two given values of W, a 
reasonable value of the enhanced moment po* 


JAFFE 


should be calculable according to the relation 


wo* = (4.606DRT 7A login k)/eze cos HAW. (29) 


In the following calculations cos #) was always 
chosen as +1 (see Section IIT). 

We may test Eq. (21) for the dielectric de- 
pendence of the rate by transforming it into the 
dimensionless form given by Eq. (28), and 
plotting W’ against 1/D?. Various theoretical 
curves will be obtained depending on the value 
of the parameter ¢ chosen. The curves in Fig. 2 
are for the two values ¢=3.2 and ¢=6.0. Experi- 
mental data agreeing with theory should fit a 
theoretical curve depending on reasonable values 
of 79 and po*. 

There is not a great amount of data available 
for testing Eq. (21). It is unfortunate that no 
author has investigated thoroughly the influence 
of both the ionic atmosphere and the dielectric 
constant of the solvent upon the rate of any 
reaction between ions and dipoles. 


(a) Dependence on Concentration 


In Fig. 1 are plotted on curve 1 the data of 
Fales and Morrell® for the inversion of sucrose 
in various concentrations of HCI, and the data 
of Kautz and Robinson® for the inversion of 
sucrose in 0.114 HCl containing various concen- 
trations of different salts. In curve 2 of Fig. 1 
the part of curve 1 up to an ionic strength of 0.1 
has been reproduced and enlarged. In curve 2 
has also been plotted the data of Brénsted and 
Grove! for the hydrolysis of ethylene acetal by 
perchloric acid. The theoretical curves average 
these points very well. It might be mentioned 
that in the case of the data of Fales and Morrell 
the points corresponding to the lowest acid con- 
centrations (below about 0.01M Cnu+) were 
omitted since they were entirely inconsistent 
with the rest of the data. In the case of inver- 
sion, observed by Kautz and Robinson, in the 
presence of NaCl two points are seen to be 
evidently in error. 

In column 2 of Table I are given the 7o values 


8H. A. Fales and J. C. Morrell, J. Am. Chem. Soc. 44, 
2071 (1922). 

9C. F. Kautz and A. L. Robinson, J. Am. Chem. Soc. 
50, 1022 (1928). 

10 J, N. Brénsted and Clinton Grove, J. Am. Chem. Soc. 
52, 1397 (1930). 
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used in fitting the various data to the theoretical 
curves in Fig. 1. It will be observed that these 
values are always reasonable and are always 
identical for the same type of atmosphere. In 
column 3 of this table are given the enhanced 
moments calculated according to Eq. (29). 
Column 4 contains the value of n? (the internal 
dielectric constant) necessary to reduce the en- 
hanced moment to 3.4X10-'§ in the case of 
sucrose! and to 410!’ in the case of ethylene 
acetal. No value of the moment for ethylene 
acetal was found recorded in literature so the 
value 4X 10-'8 was arbitrarily chosen. In the fifth 
column of this table is a key as to the value of 
the ionic strength to which this theory held for 
each charge type of ionic atmosphere. It will be 
observed that for sucrose inversion, when the 
atmosphere was made up of univalent ions, 
the theory held to a value of ionic strength >3 
which is as far as the data extended. However, 
when there were bivalent ions in the atmosphere, 
the theory applied up to a value of the ionic 
strength of about 0.6. This is similar to the case 
of reactions between ions where the presence of 
higher valence ions in the ionic atmosphere 
causes a departure from theory at lower ionic 
strengths than is the case for an ionic atmosphere 
of univalent ions. The data for the hydrolysis 
of ethylene acetal extended only to an ionic 
strength of about 0.08. 


(b) Dependence on Dielectric Constant 


In two instances there have been extensive 
studies made of the influence of the dielectric 


TABLE I. Constants used in plotting the curves of Fig. 1. 








(1) (2) i (5) (6) 
2? at 
ionic 

strength 

=1 for 

curve 1 kno 


Substances ro(A) 10'8y0* = 2? 





Sucrose +HCl 
+BaCle 
Sucrose +HC1 

+CaCle 6.0 9.05 3.32 
Sucrose +HC1 


6.0 6.65 1.91 3.867 4.02 X10~4 


3.867 4.02 X10~4 


+KCl 2.5 14.0 6.24 0.670 3.95 X10~4 
Sucrose +HCl 

+NaCl 2.5 10.5 4.18 0.670 4.301074 
Sucrose (0.15 M) 

+HC 3S 10.3 4.06 1.35 2.17 X10~4 
Sucrose (0.3 M) 

+HCl 3.5 748 2.28 1.35 
Ethylene acetal 


+HCI1Os 4.0 


2.16 X10~4 


7.3 (uw =4 X10718) 1.7 4.15 X107! 





1 FE, Landt, Naturwiss. 22, 809 (1934). 
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constant of the solvent upon the reaction rate. 
Amis and Holmes” have studied the inversion 
of 0.1000. HCl catalyzed sucrose in ethyl alcohol- 
water and in dioxane-water solvents. Since the 
data for alcohol-water media were shown from 
energy considerations to be in opposition to the 
predictions of electrostatics, they would not be 
pertinent here. However, the data for dioxane- 
water should fall in line with the requirements of 
Eq. (21). Dyas and Hill’ investigated the muta- 
rotation of glucose by HCI in methanol-water 
solvents. That part of the mutarotation brought 
about by the oxonium ion is a reaction of the 
ion-dipole type and should meet the require- 
ments of our theory. The results of these authors 
are reduced to unit normality. 

In Fig. 2 we have plotted the data of Amis 
and Holmes and also those of Dyas and Hill. 
The theoretical curves for {= 3.2, corresponding 
to a value of ro of about 3.5A used in the calcula- 
tions on atmosphere effects, averages these data 
well. However, the values for pure water as the 
solvent at 21° and 31° were omitted from the 
graph since they were so inconsistent with 
the rest of the data that they must be grossly 
in error. The Dyas and Hill data fitted the same 
curve remarkably well but led to a ro value in 
this case of only slightly more than 1A. The 
authors felt that the 7» value for the glucose- 


2 E. S. Amis and F. C. Holmes, J. Am. Chem. Soc. 63, 
2231 (1941). 

13H. E. Dyas and CG. G. Hill, J. Am. Chem. Soc. 64, 236 
(1942). 





604 E. 3. 


TABLE II. Constants used in plotting the curves of Fig. 2. 











(1) (2) (3) (4) (5) (6) 
Temp. 

Substances ro(A) 10'8y0* =? ~~ Re-o 
Sucrose +HC1 3.5 14.6 6.59 21 4.95 X1074 
Sucrose +HC1 3.5 14.7 6.65 31 2.00 X10-% 
Sucrose +HCl 3.5 14.3 6.42 41 6.76 X10-3 
Glucose +HCl 2.0 8.07 2.04 (w®=2 X10718) 20 2.95 X1071 
Glucose +HC1 2.0 8.13 2.07 (wo =2 K10718) 29.2 7.59 X107 








oxonium reaction would conceivably be less than 
the ro value for the sucrose-oxonium reaction. 
However, the difference between 1.0A and 3.5A 
for ro seemed rather large, so the choice of an 79 
value of 2.0A gave the very satisfactory results 
shown by the curve for ¢=6.0 in Fig. 2. 

In Table II, column 3 are given the enhanced 
moments calculated for the sucrose and glucose 
molecules. The values of uo* for sucrose are of 
the same order of magnitude and compare 
favorably in actual value with those obtained 
from the concentration data of entirely different 
authors when the same value of 79 is used in the 








AMIS AND G. JAFFE 


two cases. The yo* values for both sugars are 
very consistent at the different temperatures 
which is in accordance with theory because in 
the approximation (n?/2D)<1 used by us the 
ratio po*/uo depends only on x? and, therefore, 
ought to be independent of temperature. The 
values of uo* in Table I are less consistent ; these 
variations may be due to a change in the nature 
of the molecule in different surroundings. 

Column 4 of Table II contains the value of 7? 
for the moment of sucrose to be 3.4 107!§ and 
that of glucose 2X10~'§. The latter value was 
chosen arbitrarily as no value for the moment of 
glucose was found recorded in the literature. 
The values of n? are very consistent for either 
molecule, and those values for sucrose compare 
favorably with values found by the rate de- 
pendence upon ionic strength. 

It might be pointed out that our derivation of 
the potential y,; in Eq. (8) contains a theory for 
the activity coefficient of molecular dipoles, 
analogous to that of Debye and Hiickel for ions. 
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Internal Rotation in Molecules with Two 
or More Methyl Groups 
KENNETH S. PITZER 
Department of Chemistry, University of California, Berkeley, California 
July 31, 1942 

N treating the internal rotation of molecules with two 

or more identical rotating groups, it has been com- 
monly assumed, of necessity, that the potential energy for 
one rotation is independent of the position of the other 
groups. For some time the writer has been attempting to 
study this assumption but until recently the results have 
been disappointingly negative in character. However, one 
positive result may now be presented together with a brief 
discussion of the other work. 

Kennedy, Sagenkahn, and Aston! have recently pub- 
lished measurements leading to the entropy of dimethyl 
ether. After subtracting contributions for the other modes 
of motion, the value 2.65 cal. per degree at 200°K was 
obtained for the entropy contribution of the two internal 
rotational degrees of freedom. This was interpreted on the 
basis of independent potential barriers for each methyl 
group of magnitude 3100 cal. per mole. However, Anan- 
thakrishnan,? has reported Raman shifts of 160 cm™ and 
300 cm™ for dimethyl ether. Application of the usual 
harmonic oscillator formula to these frequencies yields an 
entropy contribution of 2.63 cal. per degree at 200°K in 
essentially perfect agreement with the above calorimetric 
value. Furthermore, these Raman lines have never been 
explained on any other basis. Therefore it seems likely 
that these lines are the fundamental frequencies of the 
two modes of torsional oscillation. 

However, this picture is by no means certain. Earlier 
investigators® failed to find one or both of the lines. Al- 
though the selection rules permit these torsional modes in 
the Raman effect, they predict depolarized lines whereas 
Ananthakrishnan reports both lines polarized. This latter 
difficulty is partially offset by the fact that depolarized 
lines are normally broad and both of these lines are re- 
ported as broad. Thus additional study is needed before 
this assignment can be regarded as proved or disproved. 

So far as the writer is aware, if this assignment proves 
to be correct it will be the first set of spectroscopic fre- 
quencies for torsional oscillations which is consistent with 
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reliable thermodynamic data. Assuming moderate anhar- 
monicity, these frequencies are in agreement with the 
observed entropy at higher temperatures and with the 
observed heat capacity in the gaseous state.‘ 

Figure 1 shows first the expected equilibrium configura- 
tion and then moderate displacements in each of the two 
normal modes of internal torsion. It is quite apparent that 
mode A rapidly brings hydrogen atoms very close together 
while mode B allows them to remain at considerable dis- 
tance from one another. Thus it is not surprising to find 
the first frequency almost double the second. A normal 
coordinate calculation yields for the force constants: 
2V/d¢,2=8.0X10-", 0V/dgi0g2= —3.4K10% erg per 
radian’. This value of the first constant corresponds to a 
cosine potential barrier of about 2600 cal. per mole. It is 
quite evident that the cross derivative or interaction term, 
although somewhat smaller, is by no means zero. 

The other two types of calculations concerning the inter- 
actions of two or more methyl groups will now be described 
briefly. The first concerns the possibility of using the 
entropy and heat capacity of a single gas to determine both 
the potential barrier for the rotation of one methyl group 
by itself and the magnitude of the interaction. Calculations 
for several typical cases showed that there was practically 
no chance of obtaining information by this method. The 
differences were too small to be detected experimentally. 
This means, inversely, that when potential barriers ob- 
tained from entropies are used only to calculate heat 
capacities, no serious error is likely to arise from the 
neglect of cross terms. 

The second calculation concerned the series propane, 
isobutane, and tetramethylmethane. A basic entropy was 
calculated for each molecule using the potential barrier 
from ethane (2800 cal. per mole) and the best available 
values for vibration frequencies, moments of inertia, etc. 
The amounts by which the observed entropies were smaller 
than these basic values might then be ascribed either to 
(1) interaction terms such as appear to be present in 
dimethyl! ether, or to (2) a rise in the separate potential 
barriers upon substitution of methyl groups for the hydro- 
gens in ethane. However, on the basis of reasonable ap- 
proximations, both assumptions led to the same relative 
values, 1 to 3 to 6, for the entropy differences in propane, 
isobutane, and tetramethylmethane.’ These ratios agree 
with the observed values within their respective experi- 
mental errors. However, this tells us nothing about the 
presence or absence of the interaction terms. 





A B 


Fic. 1. Left, the equilibrium position; center, moderate displacement 
in torsional mode A; right, moderate displacement in mode B. 
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We may then draw the tentative conclusion that thermo- 
dynamic data alone will not give definite information con- 
cerning cross terms in the internal rotation potential 
energy, but that spectroscopic data, particularly when 
coupled with thermodynamic data, have interesting possi- 
bilities. 


1R. M. Kennedy, M. Sagenkahn, and J. G. Aston, J. Am. Chem. 
Soc. 63, 2267 (1941). 

2R. Ananthakrishnan, Proc. Ind. Acad. Sci. A5, 285 (1937). 

3S. C. Sikar, Ind. J. Phys. 7, 257 (1932); A. Dadieu and kK. W. F. 
Kohlrausch, Monatsh. f. Chemie 57, 225 (1931); K. W. F. Kohlrausch, 
neat 68, 349 (1936); M. Wolkenstein and J. K. Syrkin, Nature 139, 288 

). 

4G. B. Kistiakowsky and W. W. Rice, J. Chem. Phys. 8, 618 (1940). 

5 This calculation may be explained very roughly as follows: (1) The 
number of interaction terms will be the number of edges of the tetra- 
hedron about the center carbon atom which connect two methyl groups. 
(2) The increments of potential barrier rise will be, respectively, 1, 2, 
and 3 for the series, but the number of degrees of freedom are, respec- 
tively, 2, 3, and 4. The products of these numbers are thus just double 
the numbers quoted in the text. 





Force Constants and Internuclear Distance 


GEORGE GLOCKLER AND GEORGE E. EVANS 
University of Iowa, Iowa City, Iowa 
July 20, 1942 


another connection we desired to know the relation 
between the force constant of the hydrogen-fluorine 
bond and its internuclear distance not only in HF but 
also in some of its more complex compounds. An attempt 
to apply Badger’s rule! indicated that the modern data of 
band spectroscopy as summarized by Herzberg? permit a 
very satisfactory estimate to be made of the internuclear 
distance if the force constant is known. The third-power 
rule was applied to the hydrides of the elements of the 
first row of the periodic table as shown in Fig. 1 and Table I. 
The results are accurate to better than one percent. The 
molecules studied were divided into three groups: (1) He, 
HD, Ds, and Het; (2) LiH, LiD, BeH*+, BeD*, BH, and 


TABLE I. Force constant K =A(R-e—B)~ megadyne/cm; R- in 
angstroms. 








Force constant (megadynes/cm) 
rom 





spectro- From % 
Molecule State A B scopy equation error 
He 1y,* 0.1152 0.1555 0.5735 0.5727 —0.14 
HD 2 28 (0.5728) 0.5728 
D2 bs m 0.5745 0.5718 —0.47 
Het 229° 1 os (0.1559) 0.1557 
Li?H 1y+ 0.1187 0.5443 (0.1021) 0.1021 
Li7D “is “ iy 0.1022 0.1023 —0.10 
(Be®H)* 4 i “ (0.2624) 0.2620 
(Be®D)* 2 = si 0.2621 0.2629 —0.31 
BUH = on i 0.3031 0.3639 
BUD a = ° 0.3164 0.3667 
Be°H sy 0.3003 0.2385 0.2253 0.2229 +1.07 
Be*D $i "= = — 0.2231 
(B"H)* 4 ee ” 0.3211 0.3225 —0.34 
CeYH Tl, st a 0.4348 0.4383 —0.80 
CRD * - 0.4347 0.4399 —1.19 
N¥4H s=- . si 0.6004 0.5876 +0.21 
O1lH 211i - “8 0.7727 0.7641 —1.11 
OVD “ os oe 0.7704 
(O16H) * 3=- - = 0.4855 0.6079 
F19H 1+ = xe 0.9626 0.9631 —0.05 
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BD; (3) BeH, BeD, BH*, CH, CD, NH, OH, OD, OH*, 
and HF. Even better results could have been obtained if 
the last one had been divided into two subgroups: BeH, 
BeD, BH*, CH, CD; and NH, OH, OD, OH*, and HF. 
The first group contains hydrogen atoms only, the second 
group is made up of molecules all of them in the '!* state, 
and the third group covers the remainder of the first-row 
atoms. 

The vibrational constant w, is not known for Be®D and 
OD, but the internuclear distances are, respectively, 
1.3427 and 0.969A. Hence it may be predicted from the 
empirical rule that they are 1527 and 2731 cm™, respec- 
tively. The values of w, for B"H and B"D are supposed 
to be 2366 and 1780 cm™, respectively, but they are un- 
certain. While they apparently fit into the third group, it 
appears plausible from Fig. 1 that these molecules should 
really belong to the second group. Since their internuclear 
distances are known accurately (1.2326 and 1.231A, respec- 
tively) their vibrational constants w, are most likely 2566 
and 1890 cm™!, as calculated from their force constants 
0.3639 and 0.3667 megadyne/cm which were obtained 
from the equations of Table I. On the other hand, the w, 
value for (B"H)*(2435.0 cm™), even though uncertain 
according to Herzberg, seems satisfactory since the mole- 
cule fits onto the straight line of the third group. But 
(O'8H)*+ does not fit the graph when the uncertain fre- 
quency is used to obtain the force constant. Hence it may 
be supposed that 3371 cm™ for w, of (O'H)* as calculated 
from our equation is more satisfactory. 

The ion (Be*H)* does not fit into group 3, nor do the 
Li?H and Li’D structures. This fact was really the occasion 
to separate the hydrides into several groups. It appears 
that the more detailed analysis of Badger’s rule has merit 
and that it can be used to obtain relations between force 
constant and internuclear distance in the case of first-row 
hydrides with satisfactory accuracy. It will of course be 
a matter of considerable interest to see how the predictions 
made will be fulfilled. 

1R. M. Badger, J. Chem. Phys. 2, 128 (1934) and 3, 710 (1935). 


2G. Herzberg, Molecular Spectra and Molecular Structure (Prentice- 
Hall, Inc., New York, 1928), Table 36. 





